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1. Introduction 



1.1. Motivation and outline 

These lectures are devoted to the Keldysh formalism for the treatment of out- 
of-equilibrium many-body systems. The name of the technique takes its origin 
from the 1964 paper of L. V. Keldysh [1]. Among the earlier approaches that are 
closely related to the Keldysh technique, one should mention J. Schwinger [2] 
and R. P. Feynman and F. L. Vernon [3]. The classical counterparts of the Keldysh 
technique are extremely useful and interesting on their own. Among them the 
Wyld diagrammatic technique [4] and the Martin-Siggia-Rose method [5] for 
stochastic systems. 

There are a number of pedagogical presentations of the method [6-8]. The 
emphasis of these notes is on the functional integration approach. It makes the 
structure of the theory clearer and more transparent. The notes also cover modern 
applications such as the Usadel equation and the nonlinear er-model. Great atten- 
tion is devoted to exposing connections to other techniques such as the equilib- 
rium Matsubara method and the classical Langevin and Fokker-Planck equations, 
as well as the Martin-Siggia-Rose technique. 

The Keldysh formulation of the many-body theory is useful for the following 
tasks: 

• Treatment of systems that are not in thermal equilibrium (either due to the 
presence of external fields, or in a transient regime) [1, 6, 8]. 

• Calculation of the full counting statistics of a quantum mechanical observable 
(as opposed to an average value or correlators) [9, 10]. 

• As an alternative to the replica and the supersymmetry methods in the theory 
of disordered and glassy systems [11-15]. 

• Treatment of equilibrium problems, in which the Matsubara analytical con- 
tinuation may prove to be cumbersome. 

The outline of these lectures is as follows. The technique is introduced and ex- 
plained for the simplest possible system, that of a single bosonic state (harmonic 
oscillator), which is later generalized to real (phonons), or complex (atoms) 
bosonic fields. Their action and Green functions are introduced in Chapter|2] Bo- 
son interactions, the diagrammatic technique and the quantum kinetic equation 
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are treated in Chapter|3] Chapter|4]is devoted to a bosonic particle in contact with 
a dissipative environment (bath). This example is used to establish connections 
with the classical methods (Langevin and Fokker-Planck) as well with the quan- 
tum equilibrium technique (Matsubara). Fermions and fermion-boson systems 
are treated in Chapter|5] Covered topics include the random phase approxima- 
tion and the quantum kinetic equation. Non-interacting Fermions in the presence 
of quenched disorder are treated in Chapter|6]with the help of the Keldysh non- 
linear er-model. 

1.2. Closed time contour 

The standard construction of the zero temperature (or equilibrium) many-body 
theory (see e.g. [7, 16]) involves the adiabatic switching "on" of interactions at 
a distant past, and "off" at a distant future. A typical correlation function has 
the form of a time ordered product of operators in the Heisenberg representa- 
tion: C(t,t') = (0\TA(t)B(t')\0), where |0) is the ground-state (or thermal 
equilibrium state) of the interacting Hamiltonian, H. This state is supposed 
to be given by |0) = S(0, — oo)|)o, where |)o is the (known) ground-state of 
the non-interacting Hamiltonian, Hq, at t — — oo. The 5-matrix operator 
S(t,t') = e lH ot e -tH(t-t ) e -iH t describes the evolution due to the interac- 
tion Hamiltonian, H — Hq, and is thus responsible for the adiabatic switching 
"on" of the interactions. An operator in the Heisenberg representation is given 



by A(t) = [S(t,0)]*A(t)S(t,0) = S(0,t)A(t)S(t,0), where A(t) is the oper- 



ator in the interaction representation. As a result, the correlation function takes 
the form: 



C{t,t') = (| T5(-oo,0)S(0, t)A{t)S{t,t')B{t')S{t' , 0)5(0, -oo)|) 
{\TA{t)B(t')S{oo,-oo)\) a 



where one employed: q(| S(— oo, 0) = e~ lL o(| S(oo, — oo)S(— oo, 0), and in- 
terchanged the order of operators, which is always allowed under the T-operation 
(time ordering). The idea is that, starting at t = — oo at the ground (or equi- 
librium) state, |)o, of the non-interacting system and then adiabatically switch- 
ing interactions "on" and "off", one arrives at t — +oo at the state |oo). The 
crucial assumption is that this state is unique, independent of the details of the 
switching procedure and is again the ground-state, up to a possible phase factor: 



Clearly this is not the case out of equilibrium. Starting from some arbitrary 
non-equilibrium state and then switching interactions "on" and "off", the system 



o(| ,S(oo,-oo)|)o 



(1.1) 




o(| £(00,-00)1)0. 
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Fig. 1. The closed time contour C. Dots on the forward and the backward branches of the contour 
denote discrete time points. 

evolves to some unpredictable state. The latter depends, in general, on the pecu- 
liarities of the switching procedure. The entire construction sketched above fails 
since we have no knowledge of the final state. 

One would like, thus, to build a theory that avoids references to the state at 
t = +00. Since traces are calculated, one still needs to know the final state. 
Schwinger's suggestion is to take the final state to be exactly the same as the 
initial one. The central idea is to let the quantum system evolve first in the for- 
ward direction in time and then to "unwind" its evolution backwards, playing the 
"movie" in the backward direction. One ends up, thus, with the need to con- 
struct a theory with the time evolution along the two-branch contour, C, depicted 
on Fig. 11.21 Then, no matter what the state at t = +00 is, after the backward 
evolution the system returns back to the known initial state. As a result, the uni- 
tary evolution operator, U ti t> = e~ tff ' t_ * along such a closed time contour is 
always a unit operator: 



In this construction there is no switching of interactions in the future. Both 
switchings "on" and "off" take place in the past: "on" - at the forward branch of 
the contour and "off" - at the backward one. This way the absence of information 
about the t = +00 state is bypassed. There is a price to pay for such luxury: a 
doubling of degrees of freedom. Indeed at every moment of time one needs 
to specify a field residing on the forward branch as well as on the backward 
branch of the contour. As a result, the algebraic structure of the theory is more 
complicated. The difficulties may be minimized, however, by a proper choice of 
variables based on the internal symmetries of the theory. 



U c = l. 



(1.2) 
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2. Free boson systems 

2.1. Partition function 

Let us consider the simplest possible many-body system: bosonic particles oc- 
cupying a single quantum state with an energy uio. It is completely equivalent, of 
course, to a harmonic oscillator. The secondary quantized Hamiltonian has the 
form: 

H = luq a' a , (2.1) 

where a) and a are bosonic creation and annihilation operators with the commu- 
tation relation [a, a*] = 1. Let us define "partition function" as: 

z M (22 ) 
Z - Tr{p} • (2 ' 2) 

If one assumes that all external fields are exactly the same on the forward and 
backward branches of the contour, then Uc = 1 and therefore Z = 1. 

The initial density matrix p = p(H) is some operator-valued function of the 
Hamiltonian. To simplify the derivations one may choose it to be the equilibrium 
density matrix, po = cxp{—/3(H — pN)} = cxp{— /3(wo — p)a^ a\. Since 
arbitrary external perturbations may be switched on (and off) at a later time, the 
choice of the equilibrium initial density matrix does not prevent one from treating 
non-equilibrium dynamics. For the equilibrium initial density matrix: 

oo 

Tr{/5 } = J2 e- f3 ^°-^ n = [1 - pMP 1 , (23) 

n=0 

where p(u>o) — e~ /3 ( a, °~ M '. An important point is that, in general, Tr{/5} is an 
interaction- and disorder-independent constant. Indeed, both interactions and 
disorder are supposed to be switched on (and off) on the forward (backward) 



Reminder: the bosonic coherent state \<f>) ((<j)\), parameterized by a complex 
number <j>, is defined as a right (left) eigenstate of the annihilation (creation) 
operator: a\(f>) = </>|</>) {{<t)\<v = (4>\4>)- The matrix elements of a normally 

ordered operator, such as the Hamiltonian, take the form 
(4>\H(a^, a)\4>') — H(4>,(f>')((f>\(f>'). The overlap between two coherent states 
is [17] = exp{00'}. Since the coherent state basis is overcomplete, the 

trace of an operator, A, is calculated with the weight: 
Tr{i} = 7T- 1 JJ d(5R0) d(Sfy) e-l^l 2 (<f>\A\<f>). 
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parts of the contour sometime after (before) t — — oo. This constant is, therefore, 
frequently omitted - it never causes a confusion. 

The next step is to divide the C contour into [2N — 2) time steps of length St, 
such that t\ = t2N = —oo and tjy = izv+i = +°° as shown in Fig. II. 21 One 
then inserts the resolution of unity in the coherent state overcomplete basis [17] 

1= /;f^) rl ^ ww (24) 

at each point j = 1,2, ... , 2N along the contour. For example, for N = 3 one 
obtains the following sequence in the expression for Tr{Ucf>o} (read from right 
to left): 

(06 | t>_,5 t 1 5 > (05 1 1>-5, | ^> 4 > (^> 4 1 1 1 3 ) <^3 I t>+d t |^»2> (02 | ^+<5 t 1 01 ) < 01 1 /5 1 6 ) , (2.5) 

where U±s t is the evolution operator during the time interval 5 t in the positive 
(negative) time direction. Its matrix elements are given by: 

(0 J+1 |f7 ± ,J0,) ee (^|e =Fi * ( ° t,o)4 %) « (0 J+1 |0,) e =Fi*<W*)* , (2.6) 

where the last approximate equality is valid up to the linear order in 5 t . Ob- 
viously this result is not restricted to the toy example, Eq i2.l\ . but holds for 
any normally-ordered Hamiltonian. Notice that there is no evolution operator 
inserted between £jy and tjv+i- Indeed, these two points are physically indistin- 
guishable and thus the system does not evolve during this time interval. 



Exercise: show that {rf>\e Ka a \(f>') = exp |0</>'e" ; | . Putting k = — /3(u>o — /it), one finds 
(</>i|po|02Jv) = exp {<^i</> 2 jvp(wo)}. 



Combining all such matrix elements along the contour together with the expo- 
nential factors from the resolutions of unity, Eq. (I2.4> . one finds for the partition 
function i2.2\ : 

2N 

i rr 2N 



e Jj,=1 



i ■ 



(2.7) 
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where the 2N x 2N matrix «G„„, stands for: 



-1 




l-h -1 




1-h -1 




1 


-l 




l+h -1 




l+h -1 . 



(2.8) 



and h = ioj^St- It is straightforward to evaluate the determinant of such a matrix 

detfzG- 1 ] = l-pMil-h 2 ^- 1 » l-pMe^ 2 ^ 1 ' - l-p(wo),(2.9) 

where one has used that <5| A" — > if N — > oo (indeed, the assumption was 
5tN — * const). Employing Eqs. jA. U and ( 12.31 1. one finds: 



det 



Tr{po} 



1 . 



(2.10) 



as it should be, of course. Notice, that keeping the upper-right element of the 
discrete matrix, Eq. ( 12. 8L is crucial to maintain this normalization identity. 

One may now take the limit A" — > oo and formally write the partition function 
in the continuous notations, <pj — > 0(t): 



Z 



/X>00exp|i J[^{t)G- x ^{t)]dt) . (2.1!) 
where according to Eqs. ( 12. 7t and d2.8i the action is given by 



where St, 



2/V 

E 

J=2 



i( t>3 



Pj-l 



5U 



u (pj<pj-i 



St; +i 



i (d)x—p{ 



^o)<p2N ) , (2.12) 



u 



G- 1 = idt 



tj-i 



u> ■ 



±St. Thus a continuous form of the operator G 1 is: 

(2.13) 



It is important to remember that this continuous notation is only an abbreviation 
that represents the large discrete matrix, Eq. ( 12.81 . In particular, the upper-right 
element of the matrix (the last term in Eq. (I2.12» . that contains the information 
about the distribution function, is seemingly absent in the continuous notations. 

To avoid integration along the closed time contour, it is convenient to split 
the bosonic field cf>(t) into the two components 4>+{t) and 4>-{t) that reside on 
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the forward and the backward parts of the time contour correspondingly. The 
continuous action may then be rewritten as 

oc oo 

S= [(tt$ + (t)[idt-wo]<t>+(t)- [dt$-(t)[idt-wo]4>-(t), (2.14) 



where the relative minus sign comes from the reverse direction of the time inte- 
gration on the backward part of the contour. Once again, the continuous notations 
are somewhat misleading. Indeed, they create an undue impression that the <f> + (t) 
and (j>_ (t) fields are completely independent from each other. In fact, they are 
connected due to the presence of the non-zero off-diagonal blocks in the discrete 
matrix, Eq. d2.8i . 



2.2. Green functions 

One would like to define the Green functions as: 



G(t, t') = -i Vcjxj) e 45 ^ 1 <j)(t)(t>{t') ee -i{<j>{t)4>{t')) , (2.15) 



where both time arguments reside somewhere on the Keldysh contour. Notice, 
the absence of the factor Z^ 1 in comparison with the analogous definition in the 
equilibrium theory [17]. Indeed, in the present construction Z = 1. This seem- 
ingly minor difference turns out to be the major issue in the theory of disordered 
systems, Chapter|6] 

According to the general property of Gaussian integrals (see |Appendix A) , the 
Green function is the inverse of the correlator matrix G _1 , Eg. ( 12. 8L standing in 
the quadratic action. Thus, one faces the unpleasant task of inverting the large 
2A*" x 2N matrix, Eq. (12.811 . It may seem more attractive to invert the differential 
operator, Eq. (12.131 1. Such an inversion, however, is undefined due to the presence 
of the zero mode (~ e~ tuJot ). The necessary regularization is provided by the 
off-diagonal blocks of the discrete matrix. The goal is to develop a formalism 
that avoids dealing with the large discrete matrices and refers to the continuous 
notations only. 

The easiest way to proceed is to recall [17] that the Green functions are traces 
of time-ordered products of the field operators (in the Heisenberg representa- 
tion), where the ordering is done along the contour C. Recalling also that the 
time arguments on the backward branch are always after those on the forward, 
one finds: 

U T<U t<\ Tr{ t(Q (f)M Trle^'qte^^'W-^po} 
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-iwo(t-i') °° 



Tr{/3 } Tr{p } 
^r^-yl ]T(™ + l)[p(^o)] m =(n + i) c -*-o(*-*') ; 



e iu (t'-t) J2. 



m=0 

- ,-<rn ^_Tr{T[a(t)aHt')}p } 



(t)<f >+ (t')) = iG 1 (t,t') = ^ (2.16) 

= 9(t- t')iG > (t, t') + 6{t' - t)iG < (t, t') ; 



, T . ^_Tr{T[a(t)aHt')]p } 



Tr{/5 } 

= 0(t' - t)iG > {t,t') + 6(t - t')iG < (t, t'); 

where the symbols T and T denote time-ordering and anti-time-ordering cor- 
respondingly. Hereafter the time arguments reside on the open time axis t G] — 
oo, oof. The Planck occupation number n stands for n.(wo) = p{u>o)/{l— p(u>o)). 

Notice the presence of non-zero off-diagonal Green functions and 
(0+ (/>_). This is seemingly inconsistent with the continuous action, Eq. J2.14> . 
This is due to the presence of the off-diagonal blocks in the discrete matrix, that 
are lost in the continuous notations. The existence of the off-diagonal Green 
functions does not contradict to continuous notations. Indeed, [idt — uio]G >,< — 
0, while [idt — ^o]G T,T = ±5(t — t'). Therefore in the obvious 2x2 matrix 
notations G _1 o G = 1, as it should be. The point is that the inverse of the 
operator [idt — wo] is not well-defined (due to the presence of the eigenmode 
(~ exp{— iu>ot}) with zero eigenvalue). A regularization must be specified and 
the off-diagonal blocks of the discrete matrix do exactly this. 

The ^-function in Eq. J2. 16b is the usual Heaviside step function. There is 
an uncertainty, however, regarding its value at coinciding time arguments. To 
resolve it, one needs to refer to the discrete representation one last time. Since 
the fields <f> always appear one time step St after the fields (j> on the Keldysh 
contour, cf. Eq. ( 12. 61 . the proper convention is: 

G T {t,t) = G f {t,t) = G < (t,t)=n. (2.17) 

Obviously not all four Green functions defined above are independent. Indeed, a 
direct inspection shows that for t ^ t'\ 



G T + G f = G> + G< 



(2.18) 
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One would like therefore to perform a linear transformation of the fields to benefit 
explicitly from this relation. This is achieved by the Keldysh rotation. 

2.3. Keldysh rotation 
Define new fields as 

<j> cl (t) = -L (0 + (t) + 0_ (t)) ; <f> q (t) = -L (0 + ( t ) _ 0_ (*)) (2.19) 

with the analogous transformation for the conjugated fields. The subscripts "cl" 
and "q " stand for the classical and the quantum components of the fields corre- 
spondingly. The rationale for these notations will become clear shortly. First, a 
simple algebraic manipulation of Eq. 12. 161 shows that 



i<Mt)&(f )) = G af3 = ( ^[*'*') GR Y ] 



(2.20) 



G R (t,t') 


- 


-G f ~ 


-G<- 




= 6(t - 


-t')(G> 


G A (t,t') 


- 


-G f - 


hG<- 


-G>) 


= 9(t' 


-t)(G< 


G K (t,t') 


- 5(°- 


+ G f - 


hG>- 




= G< 


+ G>. 



where hereafter a, (3 = (cl, q). The cancellation of the (q, q) element of this ma- 
trix is a manifestation of identity (12. 1 8i . Superscripts R, A and K stand for the 
retarded, advanced and Keldysh components of the Green function correspond- 
ingly. These three Green functions are the fundamental objects of the Keldysh 
technique. They are defined as 

G<); 

G>);(2.21) 



In the discrete representation each of these three Green functions is repre- 
sented by an N x N matrix. Since both G < and G > are, by definition, anti- 
Hermitian (cf. Eq. H.lbt ), Eq. ( 12.2 li t implies: 

G A = ^ G Rjt qK = _[ G if]t f ( 2.22) 

where the Hermitian conjugation includes complex conjugation, as well as trans- 
position of the time arguments. The retarded (advanced) Green function a is 
lower (upper) triangular matrix. Due to the algebra of triangular matrices, a prod- 
uct of any number of upper (lower) triangular matrices is again an upper (lower) 
triangular matrix. This leads to the simple rule: 

GfoGfc..oGf = G R ; 

GfoG$o...oGf = G A , (2.23) 
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G R G* G K 





Fig. 2. Graphic representation of G R , G A , and G K correspondingly. The full line represents the 
classical field component, 4> c i , while the dashed line - the quantum component, cf> q . 



where the circular multiplication signs are understood as integrations over inter- 
mediate times (discrete matrix multiplication). At coinciding time arguments, 
one finds (cf. Eqs. M.lll and ( 12. 211 ): 

G R (t,t) + G A (t,t)=0. (2.24) 

Although in the discrete representation both G R and G do contain non-zero 
(pure imaginary, due to Eqs. 12.221 . H.241 ) main diagonals (otherwise the matrix 
G is not invertible), the proper continuous convention is: 9(0) = 0. The point 
is that in any diagrammatic calculation, G R (t,t) and G A (t,t) always come in 
symmetric combinations and cancel each other due to Eq. 12.241 . It is thus a 
convenient and noncontradictory agreement to take 9(0) = 0. 

It is useful to introduce graphic representations for the three Green functions. 
To this end, let us denote the classical component of the field by a full line and 
the quantum component by a dashed line. Then the retarded Green function 
is represented by a full-arrow-dashed line, the advanced by a dashed-arrow-full 
line and the Keldysh by full-arrow-full line, see Fig. 12.31 Notice, that the dashed- 
arrow-dashed line, that would represent the (4> q 4> q ) Green function, is identically 
zero due to identity (12. 1 81 . The arrow shows the direction from <fi a towards <^g. 

For the single bosonic state (cf. Eq. J2.16» : G y — —i(n + \) e - lu) a{t-t ) an( j 
G K = -ine""' , ' t "' where n — n((jJo) — p(loq)/(1 — p(u>o)) is the Planck 
occupation number (since the system is non-interacting the initial distribution 
function does not evolve). Therefore: 

G R (t,t') = -i9(t - t') e- iw °('-*') ; 

G A (t 1 t') = i6(t' - t) e - iw °(t-t') ; (2.25) 
G K (t,t') = -i(l + 2n(uj ))e- lu)o{t - t ' ) . 

Notice that the retarded and advanced components contain information only about 
the spectrum and are independent of the occupation number, whereas the Keldysh 
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component does depend on it. Such a separation is common for systems that are 
not too far from thermal equilibrium. Fourier transforming with respect to (t — t r ) 
to the energy representation, one finds: 

G R{A \e) = (e-iu ±iO)- 1 ; (2.26) 
G K {e) = (l + 2n(w ))(-27ri)(5(e-wo) = (I + 2n(e))(-2m)5(e - u ) ■ 
Therefore for the case of thermal equilibrium, one notices therefore that 

G K (e) = coth ^ (G R (e) - G A (e)) . (2.27) 

The last equation constitutes the statement of the fluctuation-dissipation theorem 
(FDT). As it is shown below, the FDT is a general property of thermal equilibrium 
that is not restricted to the toy example, considered here. It implies the rigid 
relation between the response and correlation functions. 

In general, it is convenient to parameterize the anti-Hermitian (see Eq. \2.22V ) 
Keldysh Green function by a Hermitian matrix F = F* , as: 

G K — G R oF-Fo G A , (2.28) 

where F = F(t' , t") and the circular multiplication sign implies integration over 
the intermediate time (matrix multiplication). The Wigner transform (see below), 
/(t, e), of the matrix F is referred to as the distribution function. In thermal 
equilibrium: /(e) = coth(e/2T). 

2.4. Keldysh action and causality 

The Keldysh rotation from the (<j>+, <t>~) field components to (<f> c i, 4>q) consider- 
ably simplifies the structure of the Green functions (cf. Eqs. \2.\6\ and J2.20» . It 
is convenient, therefore, to write the action in terms of <f> c i , 4> q as well. A simple 
way of doing it is to apply the Keldysh rotation, Eq. (12. 191 . to the continuous 
action, Eq. J2.14> . written in terms of <p + , 0_ . However, as was discussed above, 
the continuous action, Eq. (I2.14I . loses the crucial information about the off- 
diagonal blocks of the discrete matrix, Eq. d2.8i . To keep this information, one 
may invert the matrix of Green functions, Eq. J2.20i . and use the result as the 
correlator in the quadratic action. The inversion is straightforward: 

G =[ g a ) =[ [G ^n [G _ 1]K ) , (2.29) 

where the three components of the inverted Green function, labelled in advance 
as A, R and K, satisfy: 

[G- 1 ]*^ [G R W}- 1 =id t -<j ±iO; (2.30) 



16 



A. Kamenev 



[G- 1 ] = -[G*]- 1 o G K o [G*]- 1 = [G*]- 1 o F - F o [G*]- 1 , 

where parameterization J2.28l > was employed in the last line. It is easy to see 
that [G fl ] _1 and [G^]" 1 are lower and upper triangular matrices correspond- 
ingly, thus justifying their superscripts. The continuous notations may create an 
impression that [G" 1 ] = (2iO)F and thus may be omitted. One should remem- 
ber, however, that this component is non-zero in the discrete form and therefore 
it is important to acknowledge its existence (even if it is not written explicitly). 

Once the correlator Eqs. H.291 . J2.30l > is established, one may immediately 
write down the corresponding action: 

oo A 1 

— oo 1 

where it is acknowledged that the correlator is, in general, a non-local function 
of time. The Green functions, Eq. J2.20I . follow from the Gaussian integral 
with this action, by construction. Notice that the presence of [G _1 ] K = (2iO)F 
(with a positive imaginary part) is absolutely necessary for the convergence of 
the corresponding functional integral. 

The structure of the Gaussian action given by Eq. ( I2.31l l is very general and 
encodes regularization of the functional integral. Since the Keldysh component 
carries the information about the density matrix, there is no further need to recall 
the discrete representation. The main features of this structure are: 

• The cl — cl component is zero. 

This zero may be traced back to identity 12. 181 . It has, however, a much simpler 
interpretation. It reflects the fact that for a pure classical field configuration (<p q = 
0) the action is zero. Indeed, in this case cf) + = 0_ and the action on the forward 
part of the contour is cancelled exactly by that on the backward part. The very 
general statement is, therefore, that 

S[(f>ci,<t> 9 = 0]=0. (2.32) 
Obviously Eq. 12.321 is not restricted to a Gaussian action. 

• The cl — q and q — cl components are mutually Hermitian conjugated upper 
and lower (advanced and retarded) triangular matrices in the time representation. 
This property is responsible for the causality of the response functions as well 
as for protecting the cl — cl component from a perturbative renormalization (see 
below). 

• The q — q component is an anti-Hermitian matrix (cf. Eq. 12.221 ) with a 
positive-definite imaginary spectrum. It is responsible for the convergence of the 
functional integral. It also keeps the information about the distribution function. 
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As was already mentioned, these three items are generic and reproduce them- 
selves in every order of perturbation theory. For the lack of a better terminology, 
we'll refer to them as the causality structure. 

2.5. Free bosonic fields 

It is a straightforward matter to generalize the entire construction to bosonic sys- 
tems with more than one state. Suppose the states are labelled by an index k, 
that may be, e.g., a momentum vector. Their energies are given by a function 
u>k, for example cok = k 2 /(2m), where m is the mass of the bosonic atoms. 
One introduces then a doublet of complex fields (classical and quantum) for 
every state k : (4>ci(t; k), <j> q (t; k)) and writes down the action in the form of 
Eq. ( 12. 311 including a summation over the index k. Away from equilibrium, the 
Keldysh component may be non-diagonal in the index k: F = F(t, t'; k, k'). 
The retarded (advanced) component, on the other hand, has the simple form 
[G^ A )]- 1 =id t -u k . 

If k is momentum, it is more instructive to Fourier transform to real space 
and to deal with (4> c i{t\ r), 4> q {t; r)). Introducing a combined time-space index 
x = (t; r), one may write down for the action of the free complex bosonic field 
(atoms): 

S Q = JJdxdx' {4>clA q ) x ( [ G H]-1 [g-l]K ) ( 0, ) (233) 

where in the continuous notations 

[G R ^]-\x,x') = 5(x-x') (id t , + , (2.34) 

\ 2m J 

while in the discrete form it is a lower (upper) triangular matrix in time (not 
in space). The [G _1 ] K component for the free field is only the regularization 
factor, originating from the (time) boundary terms. It is, in general, non-local in 
x and x', however, being a pure boundary term it is frequently omitted. It is kept 
here as a reminder that the inversion, G, of the correlator matrix must posses the 
causality structure, Eq. j2.20t . 

In an analogous way, the action of free real bosons (phonons) is (cf. Eq. jB. 9I >): 

So = ll dxdx ' ( [D-r L, ( It L ■ (235) 

where 



[D*W]-\x, x') = 6(x -x'){- d$ + vl V r 2 ,) , 



(2.36) 
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in the continuous notations. In the discrete representations [D R ( A ^]~~ 1 are again 
the lower (upper) triangular matrices. Here too the Keldysh component, [D^ 1 ] h , 
here too is just a regularization, originating from the (time) boundary terms. It 
is kept in Eq. J2.35l > to emphasize the casuality structure of the real boson Green 
function D(x, x'), analogous to Eq. J2.20I : 

v{x,x)-y D A o j, d k = d r of-fod a , 

where F — F(t, t'\ r, r') is a symmetric distribution function matrix. 



3. Collisions and kinetic equation 

3.1. Interactions 

The short range two-body collisions of bosonic atoms are described by the local 
"four-boson" Hamiltonian: H; mt = A J2 r aJaJorOr, where index r "numerates" 
spatial locations. The interaction constant, A, is related to a commonly used s- 
wave scattering length, a s , as A = 47ra s /m [18]. The corresponding term in the 
continuous Keldysh action takes the form: 



S int = -A dr dt(H>) 2 = -A dr dt [(<M+) 2 - (<M-) 2 ] • (3-D 



It is important to remember that there are no interactions in the distant past, 
t = — oo (while they are present in the future, t = +oo). The interactions 
are supposed to be adiabatically switched on and off on the forward and back- 
ward branches correspondingly. That guarantees that the off-diagonal blocks 
of the matrix, Eq. (12.81 . remain intact. Interactions modify only those matrix 
elements of the evolution operator, Eq. ( 12.6ft . that are away from t = — oo. 
It is also worth remembering that in the discrete time form the <f> fields are 
taken one time step S t after the <f> fields along the Keldysh contour C. There- 
fore the two terms on the r.h.s. of the last equation should be understood as 
(4>+(t + St)(f>+(t)) and (<j )_(t)4>_(t + St)) correspondingly. Performing the 
Keldysh rotation, Eq. ( 12.191 1. one finds 



SmMciAq] = -A dt [<j> q <f>cl{<f> 2 ci +4> 2 q ) + c.c] , (3.2) 
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Fig. 3. Graphic representation of the two interaction vertexes of the </>| 4 theory. There are also two 
conjugated vertexes with a reversed direction of all arrows. 



where c.c. stands for the complex conjugate of the first term. The collision action, 
Eq. d3.2l >. obviously satisfies the causality condition, Eq. \232\ . Diagrammati- 
cally the action J3.2I > generates two types of vertexes depicted in Fig. l3.lK as well 
as two complex conjugated vertexes, obtained by reversing the direction of the 
arrows): one with three classical fields (full lines) and one quantum field (dashed 
line) and the other with one classical field and three quantum fields. 

Let us demonstrate that adding the collision term to the action does not vi- 
olate the fundamental normalization, Z = 1. To this end one may expand 
e iS int jjj powers of A and then average term by term with the Gaussian action, 
Eq. J2.33t . To show that the normalization, Z = 1, is not altered by the colli- 
sions, one needs to show that (Si nt ) = {S 2 nt ) = . . . = 0. Applying the Wick 
theorem, one finds for the terms that are linear order in A: (4> q (i) c i(j) 2 cl + c.c.) ~ 
[G R (t, t) + G A (t, £)] G K (t, t) = 0, and (4> q 4> c i(f> 2 + c.c) = 0. The first term van- 
ishes due to identity (12.24b . while the second one vanishes because (4> q <f> q ) = 0. 
There are two families of terms that are second order in A. The first one is 
(MdfatfMfa) 2 ) ~ G R (t',t)G A (t',t)[G K (t,t')} 2 , while the second is 
(Mcifa^'M) 2 ) ~ [G R (t,t')] 2 G R (t',t)G A (t',t), where & = cj> a {t'). 
Both of these terms are zero, because G R (t', t) ~ B(t' - t), while G A (t', t) — 
G R (t,t')* ~ 6(t — t') and thus their product has no support It is easy to 
see that, for exactly the same reasons, all higher order terms vanish and thus the 
normalization is unmodified (at least in a perturbative expansion). 



'Strictly speaking, G R (t',t) and G A (t',t) are both simultaneously non-zero at the diagonal: 
t = t' . The contribution of the diagonal to the integrals is, however, ~ S?N — » 0, when N —* oo. 
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Fig. 4. Graphic representation of the two interaction vertexes of the ip 3 theory. Notice the relative 
factor of one third between them. 



As another example, consider the real boson field, Eq. J2.35i . with the cubic 
nonlinearity: 

oo oo 

S mt = ^JdrJdti P :i = ^JdrJ dt[<p 3 + -(p 3 _]=Kjdrj dt[<f&<p q + ^<pl] -(3.3) 



The causality condition, Eq. (I2.32> . is satisfied again. Diagrammatically the cu- 
bic nonlinearity generates two types of vertexes, Fig. 13. II one with two classical 
fields (full lines) and one quantum field (dashed line), and the other with three 
quantum fields. The former vortex carries the factor k, while the latter has a 
weight of k/3. Notice that for a real field the direction of the lines is not speci- 
fied by arrows. 



Exercise: Show that there are no corrections of second order in k to the partition function, Z = 1. 
Check, that the same is true for the higher orders, as well. 



3.2. Saddle point equations 

Before developing the perturbation theory further, one has to discuss the saddle 
points of the action. According to Eq. ( 12.321 , there are no terms in the action 
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that have zero power of both <f> q and cj> q . The same is obviously true regarding 
5S/5cf> c i and therefore one of the saddle point equations: 

if = ° (34) 

0<Pcl 

may always be solved by 

0,=O, (3.5) 

irrespectively of what the classical component, C /, is. One may check that this 
is indeed the case for the action given by Eqs. J2.33i plus (13.21 . Under condi- 
tion d3.5i the second saddle point equation takes the form: 

™- = ( [G R ] - 1 - A |0 C/ | 2 ) 4>d = {idt + ^ - A |0 ci | 2 ) <fe = , (3.6) 

This is the non-linear time-dependent (Gross-Pitaevskii) equation [18], that 
uniquely determines the classical field configuration, provided some initial and 
boundary conditions are specified. 

The message is that among the possible solutions of the saddle-point equa- 
tions for the Keldysh action, there is always one with a zero quantum component 
and with a classical component that obeys the classical (non-linear) equations of 
motion. We shall call such a saddle point - "classical". Thanks to Eqs. (I2.32> 
and i3.5\ . the action at the classical saddle-point field configurations is identi- 
cally zero. As was argued above, the perturbative expansion in small fluctuations 
around the classical saddle point leads to a properly normalized partition func- 
tion, Z = 1. This seemingly excludes the possibility of having any other saddle 
points. Yet, this conclusion is premature. The system may posses "non-classical" 
saddle points - such that <p q ^ 0. Such saddle points do not contribute to the par- 
tition function (and thus do not alter the fundamental normalization, Z = 1), 
however, they may contribute to the correlation functions. In general, the action 
at a non-classical saddle point is non-zero. Its contribution is thus associated 
with exponentially small (or oscillatory) terms. Examples include: tunnelling, 
thermal activation (considered in the next chapter), Wigner-Dyson level statis- 
tics, etc. 

Let us develop now a systematic perturbative expansion in deviations from the 
classical saddle point. As was discussed above, it does not bring any new infor- 
mation about the partition function. It does, however, provide information about 
the Green functions (and thus various observables). Most notably, it generates the 
kinetic equation for the distribution function. To simplify further consideration, 
let us assume that 4> c i — is the proper solution of the classical saddle-point 
equation (13.61 (i.e. there is no Bose condensate). 
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3.3. Dyson equation 

The goal is to calculate the dressed Green function, defined as: 

G al3 (t,t') = -ifvH> e l( - So+s ^ Mt)Mt') » (3-7) 



where a,(3= (cl,q) and the action is given by Eqs. d2.33i and (13.21 (or for 
real bosons: Eqs. (I2.35> and i3.3\ . with <fi — > ip). To this end one may expand 
the exponent in deviations from the classical saddle point: cf> q = and (in the 
simplest case) <j) c i — 0. The functional integration with the remaining Gaussian 
action is then performed using the Wick theorem. This leads to the standard 
diagrammatic series. Combining all one-particle irreducible diagrams into the 
self-energy matrix S, one obtains: 

G = G + Go£oG + Go±oGo±oG + ... = Go(i + £oG^, (3.8) 

where G is given by Eq. J2.20l > and the circular multiplication sign implies inte- 
grations over intermediate times and coordinates as well as a 2 x 2 matrix mul- 
tiplication. The only difference compared with the text-book [17] diagrammatic 
expansion is the presence of the 2x2 Keldysh matrix structure. The fact that 
the series is arranged as a sequence of matrix products is of no surprise. Indeed, 
the Keldysh index, a = (cl,q), is just one more index in addition to time, space, 
spin, etc. Therefore, as with any other index, there is a summation (integration) 
over all of its intermediate values, hence the matrix multiplication. The concrete 
form of the self-energy matrix, S, is of course specific to the Keldysh technique 
and is discussed below in some details. 

Multiplying both sides of Eq. (13.81 by G -1 from the left, one obtains the 
Dyson equation for the exact dressed Green function, G: 

(G- 1 -s)oG = i, (3.9) 

where 1 is the unit matrix. The very non-trivial feature of the Keldysh technique 
is that the self energy matrix, S, possesses the same causality structure as G _1 , 
Ea. (l2~29l : 



Yi A 



where are lower (upper) triangular matrices in the time domain, while 

Yi K is an anti-Hermitian matrix. This fact will be demonstrated below. Since 
both G _1 and S have the same structure, one concludes that the dressed Green 
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function, G, also possesses the causality structure, like Eq. (I2.20> . As a result, 
the Dyson equation acquires the form: 

/ [G A }-'-E A \ ( G K G R \ f 

where one took into account that [G _1 ] K is a pure regularization (~ iOF) and 
thus may be omitted in the presence of a non-zero E . Employing the specific 
form of [G fl ^] _1 , Eq. I2.34L one obtains for the retarded (advanced) compo- 
nents: 

( idt + J_ G R < A ) = E fl ( A ) o G R ( A ) . (3.12) 

\ 2m / 

Provided the self-energy component T, R ( A ^ is known (in some approximation), 
Eq. (13 . 1 21 constitutes a closed equation for the retarded (advanced) component of 
the dressed Green function. The latter carries the information about the spectrum 
of the interacting system. 

To write down the equation for the Keldysh component, it is convenient to 
parameterize it as G K = G R oF-Fo G A , where F is a Hermitian matrix in 
the time domain. The equation for the Keldysh component then takes the form: 
([G^]- 1 - E fl ) o (G R o F - F o G A ) = E g o G A . Multiplying it from the 
right by ([G" 4 ] -1 — E A ) and employing Eq. ( 13. 12L one finally finds: 

= Y, K - (Z R o F - F o E A ) , (3.13) 

where the symbol [ , ] _ stands for the commutator. This equation is the quantum 
kinetic equation for the distribution matrix F. Its l.h.s. is called the kinetic term, 
while the r.h.s. is the collision integral (up to a factor). As is shown below, 
T, K has the meaning of an "incoming" term, while T, R o F — F o T, A is an 
"outgoing" term. In equilibrium these two channels cancel each other (the kinetic 
term vanishes) and the self-energy has the same structure as the Green function: 
E A = E^oF — FoE' 4 . This is not the case, however, away from the equilibrium. 



3.4. Self-energy 

Let us demonstrate in the case of one specific example, that the self-energy ma- 
trix, E, indeed possesses the causality structure, Eq. ( 13. 10i . To this end, we con- 
sider the real boson action, Eq. J2.35l >. with the mp 3 nonlinearity, Eq. i3.3\ . and 
perform the calculations up to the second order in the parameter, n. Employing 
the two vertexes of Fig. l3~T1 one finds that: 



24 



A. Kamenev 




Fig. 5. Self-energy diagrams for the </j 3 theory. 



the cl—cl component is given by the single diagram, depicted in Fig. l3.4h . The 
corresponding analytic expression is Y, cl ~ cl (t, t') = AiK^D^t, t')D A (t, t') = 0. 
Indeed, the product D R (t, t')D A (t, t') has no support (see, however, the footnote 
in section lXTV 

the cl-q (advanced) component is given by the single diagram, Fig. 13.4b . The 
corresponding expression is: 

V A (t,t') = AiK 2 D A (t,t')D K (t,t') . (3.14) 

Since T, A (t, t') ~ D A (t, t') ~ 6{t' — t), it is, indeed, an advanced (upper trian- 
gular) matrix. There is a combinatoric factor of 4, associated with the diagram 
(4 ways of choosing external legs x 2 internal permutations x 1/(2!) for having 
two identical vertexes). 

the q-cl ( retarded) component is given by the diagram of Fig. l3.4t : 

Z R (t,t')=4iK 2 D R (t,t')D K (t,t'), (3.15) 

that could be obtained, of course, by the Hermitian conjugation of Eq. (I3.14> with 
the help of Eq. d2~22l : S fl = [Z A ] f . Since T, R (t,t') ~ D R (t,t') ~ 6{t - if), it 
is, indeed, a retarded (lower triangular) matrix. 

the q-q (Keldysh) component is given by the three diagrams, Fig. l3.4B -f. The 
corresponding expressions are: 

= 2i K 2 [D K (t,t')] 2 + 6i^y[D A (t,t')]\6iK^[D R (t,t')] 2 
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= 2iK 2 ([D K {t,t')} 2 + [D R {t,t') - D A (t,t')] 2 ^ , (3.16) 

where the combinatoric factors are: 2 for diagram d and 6 for e and f. In the last 
equality, the fact that G R (t,t')G A (t,t') = 0, due to the absence of support in 
the time domain, has been used again. Employing Eq. J2.22I . one finds Y, K — 
— [^ K ] ■ This completes the proof of the statement that £ possesses the same 
structure as D -1 , One may check that the statement holds in higher orders as 
well. In Eqs. (13.1 4i — J3 . 1 6i one has omitted the spatial coordinates, that may be 
restored in an obvious way. 



Exercise: Calculate the self-energy matrix for the </i| 4 theory to the second order in A. Show that it 
possesses the causality structure. 



3.5. Kinetic term 

To make further progress in the discussion of the kinetic equation it is convenient 
to perform the Wigner transformation (WT). The WT of a distribution function 
matrix, F(t, t'; r, r'), is a function: f (r, e; p, k), where r and p are the "center of 
mass" time and coordinate correspondingly. According to definition J2.28> . the 
F matrix appears in a product with G R — G A (or D R — D A ). Since the latter is 
a sharply peaked function at e = u>k (cf. Eq. J2.26i for free particles, while for 
interacting systems this is the condition for having well-defined quasi-particles), 
one frequently writes f (t, p, k), understanding that e = lo^. 

To rewrite the kinetic term (the l.h.s. of Eq. J3.131 ) in the Wigner represen- 
tation, one notices that the WT of idt is e, while the WT of is —k 2 . Then 



Reminder: The Wigner transform of a matrix A(r, r') is defined as 




One may show that the Wigner transform of the matrix C = A o B is equal to: 

c(p, k) =JJ dr ^JJ a (p + y> k + b (p + y> k + e^~ k ^). 
Expanding the functions under the integrals in ki and r^, one finds: 
c(p, k) = a(p, k) b(p, k) + (2i)~ 1 (V p aV k b - V fe aV p fe) + . . . . 
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e.g. [F, V 2 .]- — > [/c 2 ,f]- + iVfefc 2 V p f = 2ikV p f, where the commutator van- 
ishes, since WT's commute. In a similar way: [F,i<9 t ]_ — > —id T i. If there 
is a scalar potential V^(r)aja r in the Hamiltonian, it translates into the term 
-V(4> c i4> q + 4> q (f>ci) in the action and thus — V(r) is added to [G^^] -1 . This, 
in turn, brings the term — [F, V]_ to the l.h.s. of the Dyson equation 13.131 . or 
after the WT: iEV \f, where E = —V P V is the electric field. As a result, the 
WT of the Dyson equation ( 13.91 takes the form: 

(ar-«fcVp-SV fc )f(r,p,fc)= Jeoilf], (3.17) 

where Vk = fe/m and I co ;[f] is the WT of the r.h.s. of Eq. 13. 131 (times i). This 
is the kinetic equation for the distribution function. 

For real bosons with the dispersion relation e = ujk, the kinetic term is (cf. 
Eq. OJ): [e 2 - ^fc,F]_ -» 2z(ed T - w fc (V fc w fe )V p )f = 2ie(d T - v k V 
where = V k^k is the group velocity. As a result, the kinetic equation takes 
the form: (d r — Vk V p ) f (r, p, k) — I co i [f], where the collision integral I co i [f] is 
the WT of the r.h.s. of Eq. dTlH . divided by -2ie. 

3.6. Collision integral 

Let us discuss the collision integral, using the (p 3 theory calculations of section 
I3.4l as an example. To shorten the algebra, let us consider a system that is spatially 
uniform and isotropic in momentum space. One, thus, focuses on the energy 
relaxation only. In this case the distribution function is f (t, p, k) — f (r, u>k) — 
f (t, e), where the dependence on the modulus of the momentum is substituted by 
the Uk = e argument. Employing Eqs. (13. 14-1 — J3. 161 . one finds for the WT of the 
r.h.s. of Eq. ( l3~L3l 2 : 

E K o F - F o E- 4 -> -2i f (t, e) ft du; M(r, e, w) (f (r, e - w) + f (r, w)) ; 

S K -> -2ifdu M(T,e,u)(f(T,e - u)f{T,u) + 1) , (3.18) 

where the square of the transition matrix element is given by: 

M(r, e, w) = 2ttk 2 ^ A d (r, e-uj-k-q) A d (r, 9 ) . (3.19) 

9 

2 Only products of WT's are retained, while all the gradient terms are neglected, in particular 
D K — > f (d R — d A ). The energy-momentum representation is used, instead of the time-space 
representation as in Eqs. 13. 141 - 13. 161 . and in the equation for T, R o F — F o T, A one performs a 
symmetrization between the ui and e — w arguments. 
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Here A d = i(d R - d A )/(2vr) and d R ( A )(r, e; k) is the WT of the retarded 
(advanced) Green function. One has substituted the dressed Green functions into 
Eqs. ( 13 . 1 4-1 — J3 . 1 61 instead of the bare ones to perform a partial resummation of 
the diagrammatic series. (This trick is sometimes called the self-consistent Born 
approximation. It still neglects the vertex corrections.) Assuming the existence of 
well defined quasi-particles at all times, one may regard Aa(r, e, k) as a sharply 
peaked function at e = lo^. In this case Eq. ( 13.191 1 simply reflects the fact that 
an initial particle with e = u>k decays into two real (on mass-shell) particles with 
energies uj — uj q ande — u) — u>k-q- As a result, one finally obtains for the kinetic 
equation: 



where T is a temperature. This is the thermal equilibrium distribution function. 
According to the kinetic equation (I3.20> . it is stable for any temperature (the latter 
is determined either by an external reservoir, or, for a closed system, from the 
total energy conservation). Since the equilibrium distribution obviously nullifies 
the kinetic term, according to Eq. ( I3.13> the exact self-energy satisfies T, K = 
coth(e/(2T))(S fl — T, A ). Since also the bare Green functions obey the same 
relation, Eq. ( 12.271 . one concludes that in thermal equilibrium the exact dressed 
Green function satisfies: 



This is the statement of the fluctuation-dissipation theorem (FDT). Its conse- 
quence is that in equilibrium the Keldysh component does not contain any ad- 
ditional information with respect to the retarded one. Therefore, the Keldysh 
technique may be, in principle, substituted by a more compact construction - the 
Matsubara method. The latter does not work, of course, away from equilibrium. 

Returning to the kinetic equation ( I3.20L one may identify "in" and "out" 
terms in the collision integral. Most clearly it is done by writing the collision 
integral in terms of the occupation numbers rifc, defined as f = 1 + 2n. The 
expression in the square brackets on the r.h.s. of Eq. (I3.20> takes the form: 
4 [n e _ w n w — n e (n e _ w + n w + 1)]. The first term: n c _ w n w , gives a probabil- 
ity that a particle with energy e — w absorbs a particle with energy u> to populate a 




f (e - w)f (w) + 1 - f (e) (f (e + f (w)) , (3.20) 



where the time arguments are suppressed for brevity. 

Due to the identity: coth(a- b) coth(fe) + 1 = coth(a) (coth(a- b) +coth(6)), 
the collision integral is identically nullified by 





(3.22) 
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state with energy e - this is the "in" term of the collision integral. It may be traced 
back to the Y, K part of the self-energy. The second term: — n e (n e _ w + n u + 1), 
says that a state with energy e may be depopulated either by stimulated emission 
of particles with energies e — uo and u>, or by spontaneous emission (unity). This 
is the "out" term, that may be traced back to the Y, R ( A '> contributions. 

Finally, let us discuss the approximations involved in the Wigner transforma- 
tions. Although Eq. J3 . 1 31 is formally exact, it is very difficult to extract any 
useful information from it. Therefore, passing to an approximate, but much more 
tractable, form like Eqs. (I3.17> or J3.20i is highly desirable. In doing it, one has 
to employ the approximate form of the WT. Indeed, a formally infinite series 
in VfcVp operators is truncated, usually by the first non-vanishing term. This 
is a justified procedure as long as 5k Sp ^> 1, where 5k is a characteristic mi- 
croscopic scale of the momentum dependence of f , while 5p is a characteristic 
scale of its spatial variations. One may ask if there is a similar requirement in the 
time domain: 5e 5t ^ 1, with 5e and 5t being the characteristic energy and the 
time scale of f, correspondingly? Such a requirement is very demanding, since 
typically 5e T and at low temperature it would allow to treat only very slow 
processes: with 5t ^> 1/T. Fortunately, this is not the case. Because of the 
peaked structure of Ad(e, k), the energy argument e is locked to uj^ and does 
not have its own dynamics as long as the peak is sharp. The actual criterion is 
therefore that 5e is much larger than the width of the peak in Aj(e, k). The 
latter is, by definition, the quasi-particle life-time, r qp , and therefore the condi- 
tion is T qp ^> 1/T. This condition is indeed satisfied by many systems with the 
interactions that are not too strong. 



4. Particle in contact with an environment 

4.1. Quantum dissipative action 

Consider a particle with the coordinate living in a potential [/($) and at- 
tached to a harmonic string ip(t; x). The particle may represent a collective de- 
gree of freedom, such as the phase of a Josephson junction or the charge on a 
quantum dot. On the other hand, the string serves to model a dissipative envi- 
ronment. The advantage of the one-dimensional string is that it is the simplest 
continuum system, having a constant density of states. Due to this property it 
mimics, for example, interactions with a Fermi sea. A continuous reservoir with 
a constant density of states at small energies is sometimes called an "Ohmic" en- 
vironment (or bath). The environment is supposed to be in thermal equilibrium. 
The Keldysh action of such a system is given by the three terms (cf. Eqs. ( IB. 5> 
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and <|235}): 

Sp[Q]= I dt 



-2$ 



q dt 2 



(it / <ix <£> T l) 1 ; 



(4.1) 



where we have introduced vectors of classical and quantum components, e.g. 
$ T = ($ c /,$ 9 ) and the string correlator, D~ l , is the same as in Eqs. J2.35t . 
(I2.36> . The interaction term between the particle and the string is taken to be the 
local product of the particle coordinate and the string stress at x — (so the force 
on the particle is proportional to the local stress of the string). In the time domain 
the interaction is instantaneous, ^(t)'V x (p(t,x)\ x —o — > $+Vy>+ — $_V</?_ on 
the Keldysh contour. Transforming to the classical-quantum notations leads to: 
2(<l> c ;V(^ g + &q^7ip c i), that satisfies the causality condition, Eq. (|T 



2J. In the 

matrix notations it takes the form of the last line of Eq. ( 14. U . where <j\ is the 
standard Pauli matrix. The interaction constant is ypy. 

One may now integrate out the degrees of freedom of the Gaussian string to re- 
duce the problem to the particle coordinate only. According to the standard rules 
of Gaussian integration (see. |Appendix A) , this leads to the so-called dissipative 
action for the particle: 



Sdi. 



-j J J dtdt' $ 2 (t) a{ V x V x ,D(t 



x—x'—O 



ai $(«')■ (4.2) 



i(t-t') 



The straightforward matrix multiplication shows that the dissipative correlator 
L _1 possesses the standard causality structure of the inverse Green function, e.g. 
Eq. ( I2.29> . Fourier transforming its retarded (advanced) components, one finds: 



L R(A) {e 



E 



k' 2 



{e±iO) 2 - k 2 



± — e + const , 



(4.3) 



where we put v s = 1 for brevity. The e-independent constant (same for R and 
A components) may be absorbed into the redefinition of the harmonic part of the 
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potential U ($) = const <f> 2 + . . . and, thus, may be omitted. In equilibrium the 
Keldysh component of the correlator is set by the FDT: 

[L-'] K (e) = coth^ {[L*]- 1 - [L*]- 1 ) = iecoth^ . (4.4) 

It is an anti-Hermitian operator with a positive-definite imaginary part, rendering 
convergence of the functional integral over <1>. 

In the time representation the retarded (advanced) component of the correlator 
takes a simple local form: [i fl '' 4 )] = =pi 5(t — t') d t >. On the other hand, at 
low temperatures the Keldysh component is a non-local function, that may be 
found by the inverse Fourier transform of Eq. 14.41 : 



Finally, for the Keldysh action of the particle connected to a string, one obtains: 



= Jdt 



-^i'. / (^ + ^J-^(«^ + «i',)+r(«iv/-«i',i 



fJ*"*M ^£(t-i<)) *M' (4 ' 6) 



V) 



This action satisfies all the causality criterions listed in section l2~4l Notice, that in 
the present case the Keldysh (q — q) component is not just a regularization factor, 
but rather a quantum fluctuations damping term, originating from the coupling 
to the string. The other manifestation of the string is the presence of the friction 
term, ~ jdt in the R and the A components. In equilibrium the friction coeffi- 
cient and fluctuations amplitude are rigidly connected by the FDT. The quantum 
dissipative action, Eq. ( 14. 6> . is a convenient playground to demonstrate various 
approximations and connections to other approaches. 

4.2. Saddle-point equation 

The classical saddle point equation (the one that takes <& g (t) — 0) has the form: 



1 6S[$\ 



2 5<t> n 



_ d^ cl 7 d$ d dU&d) _ 

dt 2 2 dt 9$ c/ } 

* o =0 



This is the deterministic classical equation of motion. In the present case it hap- 
pens to be the Newton equation with the viscous force: — ( r y/2)$> c i. This approx- 
imation neglects both quantum and thermal fluctuations. 
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4.3. Classical limit 

One may keep the thermal fluctuations, while completely neglecting the quantum 
ones. To this end it is convenient to restore the Planck constant in the action ( 14. 6> 
and then take the limit h —* 0. For dimensional reasons, the factor h^ 1 should 
stand in front of the action. To keep the part of the action responsible for the 
classical equation of motion J4.7i free from the Planck constant it is convenient 
to rescale the variables as: <I> g — > TiQ q . Finally, to have temperature in energy 
units, one needs to substitute T — > T/% in the last term of Eq. J4.6I >. The limit 
?i — > is now straightforward: (i) one has to expand U($ c i ± Ti$ q ) to the first 
order in Ti$ q and neglect all higher order terms; (ii) in the last term of Eq. J4.6I 
the ft — > limit is equivalent to the T — > oo limit, see Eq. J4.5I >. As a result, the 
classical limit of the dissipative action is: 

— oo 

Physically the limit h — ► means that hil -C T, where fl is a characteristic 
classical frequency of the particle. This condition is necessary for the last term 
of Eq. ( 14. 6> to take the time-local form. The condition for neglecting the higher 
order derivatives of U is fi <C J^l' where $cl is a characteristic classical ampli- 
tude of the particle motion. 



4.4. Langevin equations 

One way to proceed with the classical action J4.8i is to notice that the exponent 
of its last term (times i) may be identically rewritten in the following way: 



Jv£(t) e^fe^W^W^'W] . (4.9) 



This identity is called the Hubbard-Stratonovich transformation, while is an 
auxiliary Hubbard-Stratonovich field. The identity is proved by completing the 
square in the exponent on the r.h.s., performing the Gaussian integration at every 
instance of time and multiplying the results. There is a constant multiplicative 
factor hidden in the integration measure, X>£. 

Exchanging the order of the functional integration over £ and one finds for 
the partition function: 
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Since the last (imaginary) exponent depends only linearly on & q (t), the integra- 
tion over results in the ^-function of the expression in the round brackets. 
This functional <5-function enforces its argument to be zero at every moment of 
time. Therefore, among all the possible trajectories $ c ;(t), only those that satisfy 
the following equation contribute to the partition function: 

d 2 $ c/ 7# d dU($ c i) _ . 

This is a Newton equation with a time dependent external force £(t). Since, the 
same arguments are applicable to any correlation function of the classical fields, 
e.g. ($ci(t)Q c i (t')), a solution strategy is as follows: (i) choose some realization 
of £(t); (ii) solve Eq. (I4.1H (e.g. numerically); (iii) having its solution, $ c /(i), 
calculate the correlation function; (iv) average the result over an ensemble of 
realizations of the force £(t). The statistics of the latter are dictated by the weight 
factor in the T>£ functional integral. It states that £(t) is a Gaussian short-range 
(white) noise with the correlators: 

(£(*)> =0; (Z(t)tV))=<TT5(t-lf). (4.12) 

Equation (14.111 with the white noise on the r.h.s. is called the Langevin equation. 
It describes classical Newtonian dynamics in presence of stochastic thermal fluc- 
tuations. The fact that the noise amplitude is related to the friction coefficient, 7 
and to the temperature is a manifestation of the FDT The latter holds as long as 
the environment (string) is at thermal equilibrium. 

4.5. Martin-Siggia-Rose 



In section l4T4l one derived the Langevin equation for a classical coordinate, $ c /, 
from the action written in terms of $ c / and another field, $ 9 . An inverse pro- 
cedure of deriving the effective action from the Langevin equation is known as 
the Martin-Siggia-Rose (MSR) [5] technique. It is sketched here in the form 
suggested by De-Dominics [5]. 
Consider a Langevin equation: 

d[#]=f(t), (4.13) 

where is a (non-linear) differential operator acting on the coordinate <£>(t) 
and £(t) is a white noise force, specified by Eq. (14. 12i . Define the "partition 
function" as: 

Z[£] = Jv*J[6]5(6[$\-£tt))=\. (4.14) 
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It is identically equal to unity by virtue of the integration of the 5-function, pro- 
vided J\(D\ is the Jacobian of the operator O [$] . The way to interpret Eq. J4. 14i 
is to discretize the time axis, introducing TV-dimensional vectors <£>j = $(tj) 
and £j — £,(tj). The operator takes the form: O l — + Tijk^j^k + • • •, 

where a summation is taken over repeated indexes. The Jacobian, J, is given 
by the absolute value of the determinant of the following N x N matrix: Jy = 
dOi/d&j — Oij + 2Tijk^k + ■ ■ ■■ It is possible to choose a proper (retarded) 
regularization where the Jy matrix is a lower triangular matrix with a unity main 
diagonal (coming entirely from the On = 1 term). Clearly, in this case, J — 1. 
Indeed, consider, for example, 0[<f>] = dt& — W(G>). The retarded regularized 
version of the Langevin equation is: $^ = + 5t(W(&i-i) Clearly 
in this case Jn = 1 and Ji.i-i = — 1 — W (&i-i)5 t , while all other components 
are zero; as a result J = 1. 

Although the partition function J4. 14i is trivial, it is clear that all the meaning- 
ful observables and the correlation functions may be obtained by inserting a set of 
factors: ... in the functional integral, Eq. J4. 14i . Having this in mind, 

let us proceed with the partition function. Employing the integral representation 
of the 5-function with the help of an auxiliary field ^(t), one obtains: 




(4.15) 



where O stands for the retarded regularization of the O operator and thus one 
takes J = 1. One may average now over the white noise, Eq. (14. 12i . by per- 
forming the Gaussian integration over £: 




The exponent is exactly the classical limit of the Keldysh action, cf. Eq. (I4.8i (in- 
cluding the retarded regularization of the differential operator), where $ = $ c ; 
and $ = $ f/ . The message is that the MSR action is nothing, but the classical 
(high temperature) limit of the Keldysh action. The MSR technique provides a 
simple way to transform from a classical stochastic problem to its proper func- 
tional representation. The latter is useful for an analytical analysis. One example 
is given below. 

4.6. Thermal activation 

Consider a particle in a meta-stable potential well, plotted in Fig. 14. 6h . The 
potential has a meta-stable minimum at <f> = and a maximum at $ — 1 with the 
amplitude Uq. Let us also assume that the particle's motion is over-damped, i.e. 
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Fig. 6. a) A potential with a meta-stable minimum, b) The phase portrait of the Hamiltonian system, 
Eq. 14.191 . Thick lines correspond to zero energy, arrows indicate evolution direction. 



7 ^> \J\J" . In this case one may disregard the inertia term, leaving only viscous 
relaxation dynamics. The classical dissipative action J4.8I takes the form: 



S[$] = 2 dt 



7 dUjQd) 
q[ 2 dt d®c. 



(4.17) 



The corresponding saddle point equations are: 

2ijT<S> q ; 



2 - _ du{<s> cl ) 

2 9cl d$ cl 



7 



(4.18) 



These equations possess the classical solution: $ 9 (t) = and $ c /(t) satisfies 
the classical equation of motion: ^ $ c / = —dU($ c i)/d$> c i- For the initial con- 
dition $ c ;(0) < 1 me latter equation predicts the viscous relaxation towards the 
minimum at $ c ; = 0. According to this equation, there is no possibility to es- 
cape from this minimum. Therefore the classical solution of Eqs. ( I4.18l > does not 
describe thermal activation. Thus one has to look for another possible solution of 
Eqs. (14. 181 . the one with $ 9 ^ 0. 

To this end let us make a simple linear change of variables: $ c ;(i) = q(t) 
and 3> g (t) = p(t) I (yi). Then the dissipative action (14. 17t acquires the form of a 
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Hamiltonian action: 

9, r BTT(n\ „1 

,(4.19) 



iS = - Jdt(pq-H(p,q)); H{p,q) = - 



dU(q) _ 2 
-P — ^ 1" T P 



where the fictitious Hamiltonian, 7?, is introduced 3 . It is straightforward to 
see that in terms of the new variables the equations of motion (14. 1 81 take the 
form of the Hamilton equations: q = dH/dp and p = —dH/dq. One needs, 
thus, to investigate the Hamiltonian system with the Hamiltonian Eq. ( I4.19> . To 
visualize it, one may plot its phase portrait, consisting of lines of constant energy 
E = H(p(t),q(t)) on the (p,q) plane, Fig. 14.6b . The topology is determined 
by the two lines of zero energy: p — and Tp = dU(q)/dq, that intersect 
at the two stationary points of the potential: q = and q = 1. The p = 
line corresponds to the classical (without Langevin noise) dynamics (notice, that 
the action is identically zero for motion along this line) and thus q = is the 
stable point, while q = 1 is the unstable one. Due to Liouville theorem, every 
fixed point must have one stable and one unstable direction. Therefore, along 
the "non-classical" line: p = T^ 1 dU(q)/dq, the situation is reversed: q = is 
unstable, while q = 1 is stable. It is clear now that, to escape from the bottom 
of the potential well, q = 0, the system must move along the non-classical line 
of zero energy until it reaches the top of the barrier, q = 1, and then continue 
to drop according to the classical equation of motion (moving along the classical 
line p = 0). There is a non-zero action associated with the motion along the 
non-classical line: 



l 



iS = -Jdtpq = - jp(q)dq = ~J dq = - ^ , (4.20) 



where one has used that H = Q along the integration trajectory. As a result, the 
thermal escape probability is proportional to e lS — e~ u °/ T , which is nothing but 
the thermal activation exponent. 



4. 7. Fokker-Planck equation 

Another way to approach the action < I4. 17i is to notice that it is quadratic in Q q 
and therefore the T><fr q integration may be explicitly performed. To shorten nota- 
tions and emphasize the relation to the classical coordinate, we shall follow the 

3 Amazingly, this trick of rewriting viscous (or diffusive) dynamics as a Hamiltonian one, works 
in a wide class of problems. The price, one has to pay, is the doubling of the number of degrees 
of freedom: q and p in the Hamiltonian language, or "classical" and "quantum" components in the 
Keldysh language. 
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previous section and denote ( f ) c ;(t) = q(t). Performing the Gaussian integra- 
tion over <£> g of e lS [*l, with 5[$ c /, $ q ] given by Eq. J4.17> . one finds the action, 
depending on $ c ; = q only: 



iS[q] = 



2 7 T 



dt 



(4.21) 



One may now employ the same trick, that allows to pass from the Feynman 
path integral to the Schrtidinger equation. Namely, let us introduce the "wave 
function", V(q, t), that is a result of the functional integration of e lS ^ over all 
trajectories that at time t pass through the point q^ = q. Adding one more 
time step, S t , to the trajectory, one may write T^ft/jv, * + 8t) as an integral of 
V{q N -i,t) = V{q + 6 q ,t) overS q = q N _ 1 - q : 



V(q,t + S t ) =C [dS q e ~^(* ~k + U '« {q+Sq) Y V{q + 5 q ,t) 



(4.22) 



where the factor C from the integration measure is determined by the condition: 
Cjd5 q cxp {— 7<5^/(8T<5 t )} = 1. Expanding the expression in the square brack- 
ets on the r.h.s. of the last equation to the second order in S q and the first order in 
St, one finds: 





7 



7 



7 



(4.23) 



where (8%) = Cfd5 q exp {-j8%/(8T8 t )} 8 2 q = iT8 t /j. Finally, rewriting the 
last expression in the differential form, one obtains: 



dV 
dt 



d_dU 
dq dq 



■7; — ^ h T- 



O 2 



v-~ — 

7 9q 



dU m dV 

-p + T 

dq dq 



(4.24) 



This is the Fokker-Planck (FP) equation for the evolution of the probability dis- 
tribution function, V(q, t). The latter describes the probability to find the particle 
at the point q(= $) at time t. If one starts from an initially sharp (determin- 
istic) distribution: V(q,0) = 8(q — 9(0)), then the first term on the r.h.s. of 
the FP equation describes the viscous drift of the particle in the potential U(q). 
Indeed, in the absence of the second term (T = 0), the equation is solved by 
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V(q, t) = 6(q — q(t)), where q(t) satisfies the deterministic equation of motion 
(-f/2)q(t) — —dU(q(t))/dq 4 . The second term describes the diffusion spread- 
ing of the probability distribution due to the thermal stochastic noise £(t). For a 
confining potential U(q) (such that J7(±oo) — > oo) the stationary solution of the 
FP equation is the equilibrium Boltzmann distribution: P(q) ~ exp{— U(q)/T}. 

The FP equation may be considered as the (imaginary time) Schrodinger equa- 
tion: V = HP, where the "Hamiltonian", H, is nothing but the "quantized" ver- 
sion of the classical Hamiltonian, introduced in the previous section, Eq. J4.191 . 
The "quantization" rule is p — > p = —d/dq, so the canonical commutation re- 
lation: = 1, holds. Notice that before applying this quantization rule, 
the corresponding classical Hamiltonian must be normally ordered. Namely, 
the momentum p should be to the left of the coordinate q, cf. Eq. ( I4.19> . Us- 
ing the commutation relation, one may rewrite the quantized Hamiltonian as: 
H = Tp 2 — pU' q = T (p— U'J(2T)) (p - U'J{2T)) - {U' q f/(AT) + U'< q /2 
(we took 7/2 = 1) and perform the canonical transformation: Q = q and 
P = p — Ug/(2T). In terms of these new variables the Hamiltonian takes the 
familiar form: H = TP 2 + V(Q), where V(Q) = -{U' Q ) 2 I{AT) + U^ Q /2, 
while the "wave function" transforms as V(Q, t) = e u ^^ 2T ^P. 

4.8. From Matsubara to Keldysh 

In some applications it may be convenient to derive an action in the equilibrium 
Matsubara technique and change to the Keldysh representation at a later stage to 
tackle out-of-equilibrium problems. This section intends to illustrate how such 
transformation may be carried out. To this end consider the following bosonic 
Matsubara action: 

00 ^ 

5[* m ] = 7 r 2 i e ™n $ ™i 2 ' (4 ' 25) 



4 To check this statement one may substitute V(q, t) = 5(q — q(t)) into the T = FP equation: 
S' q (q - q(t))(-q(t)) = (2/ 7 ) [U q " q 8{q - q(t)) + U q S' q (q - q(t))] . Then multiplying both parts 
of this equation by q and integrating over dq (by performing integration by parts), one finds: q{t) = 
~{2h)U' q {q{t)). 



Reminder: The Matsubara technique deals with the imaginary time r confined to the interval 
r£ [0, f}[ , where /3 = 1/T. All bosonic fields must be periodic in this interval: </>(t + 0) = 4>(t), 
while the fermionic fields are antiperiodic: if)(r 4-/3) = — iK 1 ")- K i s convenient to introduce the 



discrete Fourier (Matsubara) transform, e.g. <j> m = J dr<j>(T) e !e m T , where for bosons 



e m = 2ivmT, while for fermions e m = ir(2m + 1)T and m = 0, ±1, . . .. 
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where $ m = <f>_ m are the Matsubara components of a real bosonic field, 3>(t). 
Notice, that due to the absolute value sign: |e TO | ^ id T . In fact, in the imaginary 
time representation the action J4.25I has the non-local form: 



Sm=-} ffdrdr'Hr) , ?J $(r') . (4.26) 

2 jj sin (ttT{t - t')) 

o 

This action is frequently named after Caldeira and Leggett [19], who used it to 
investigate the influence of dissipation on quantum tunnelling. 

To transforn to the Keldysh representation one needs to double the number 
of degrees of freedom: $ — > $ = (<E> C ;, <& g ) T . Then according to the causal- 
ity structure, section 12.41 the general form of the time translationally invariant 
Keldysh action is: 

where [L R ^ A > (e)] _1 is the analytical continuation of the Matsubara correlator 
|e m |/2 from the upper (lower) half-plane of the imaginary variable e m to the real 
axis: — ie m — > e. As a result, [L^^fe)] -1 = ±ie/2. The Keldysh component 
follows from the FDT: [L- 1 ]*^) = iecothe/(2T), cf. Eas. JO and d4T4l. 
Therefore the Keldysh counterpart of the Matsubara action, Eqs. ( 14.251 or J4.26I 
is the already familiar dissipative action, Eq. J4.6I >. (without the potential terms, 
of course). One may now include external fields and allow the system to deviate 
from the equilibrium. 

4. 9. Dissipative chains and membranes 

Instead of dealing with a single particle connected to a bath, let us now consider 
a chain or lattice of coupled particles, with each one connected to a bath. To 
this end, one (i) supplies a spatial index, r, to the field: Q(t) — > <£>(t;r), and 
(ii) adds the harmonic interaction potential between nearest neighbors particles: 
~ r) — <I>(i, r + l)) 2 — > (V r <f>) 2 in the continuous limit. By changing 
to the classical-quantum components and performing the spatial integration by 
parts (cf. Eq. IB. 9» . the gradient term translates to: <5> g V 2 $ c / + $ c / V 2 $ g . Thus 
it modifies the retarded and advanced components of the correlator, but it does 
not affect the (q — q) Keldysh component: 

[L^]- 1 = \s(t- S(r - r')( T d t , + DV 2 r ,) , (4.28) 

where D is the rigidity of the chain or the membrane. In the Fourier representa- 
tion: [L R ^ (e; k)]^ 1 = i( ± it — Dk 2 ^. In equilibrium the Keldysh component 
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is not affected by the gradient terms, and is given by Eq. J4.4I (in the real space 
representation it acquires the factor 6(r — r')). In particular, its classical limit 
is (cf. Eq. lO) [L~ X ] K = i2T5(t - t')S(r - r'). As a result, the action of a 
classical elastic chain in contact with a bath is: 



oo 

S[$] = 2 f dr f dt 



A> l (^-D^ cl + ^§^-)+r2T^ 



, (4.29) 



where the inertia terms have been neglected and we put 7/2=1 for brevity. 

One may introduce now an auxiliary Hubbard-Stratonovich field r) and 
write the Langevin equation according to section l4~4l 



$ cJ -£V^ ci + ^M=^;r), (4.30) 

where £ is a Gaussian noise: (£(t; r)£(t'; r')) = 2T5(t — t')5(r — r') with short- 
range correlations. 

Let us consider an elastic chain sitting in the bottom of the (r-independent) 
meta-stable potential well, depicted in Fig. 14.6b . If a sufficiently large piece of 
the chain thermally escapes from the well, it may find it favorable to slide down 
the potential, pulling the entire chain out of the well. To find the shape of such an 
optimally large critical domain and its action, let us change to the Hamiltonian 
variables of section l4~6l q(t; r) = $ c /(t; r) and p(t; r) = 2i$ g (i; r). The action 

.29i takes the Hamiltonian form: 



iS = -JJdrdt(pq- H(p,q)) ; H = -p + p DV 2 r q + Tp 2 , (4.31) 

and the corresponding equations of motion are: 
SH 

q = — = DV 2 r q-U' q (q) + 2Tp; (4.32) 
P =-^ = -DW 2 rP + pU^(q). 

These are complicated partial differential equations, that cannot be solved in gen- 
eral. Fortunately, the shape of the optimal critical domain can be found. As was 
discussed in section I4~6l the minimal action trajectory corresponds to a motion 
with zero energy, H = 0. According to Eq. ( 14.3 1> . this is the case if either p = 
(classical zero-action trajectory), or Tp = U' q (q) — DW 2 q (finite-action escape 
trajectory). In the latter case the equation of motion for q(t; r) takes the form of 
the classical equation in the reversed time: q = —DW 2 q + U' q (q) = Tp . Thanks 
to the last equality the equation of motion for p(t; r) is automatically satisfied 
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5 . In the reversed time dynamics the q(t; r) = configuration is unstable and 
therefore the chain develops a "tongue" that grows until it reaches the stationary 
shape: 

-DV 2 r q + U' q (q)=0. (4.33) 

The solution of this equation gives the shape of the critical domain. Once it is 
formed, it may grow further according to the classical equation q = DV 2 q — 
U' q (q) and p = with zero action. The action along the non-classical escape 
trajectory, paid to form the "tongue" is {H(p, q) = 0): 

-iS=JJdrdtpq=yjdrdt(U' q (qyDV 2 r q)q=^ Jdr(lJ(q) + j(V r q) 2 ) , (4.34) 

where in the last equality an explicit integration over time is performed. The 
escape action is given therefore by the static activation expression that includes 
both the potential and the elastic energies. The optimal domain, Eq. J4.33i . is 
found by the minimization of this static action J4.34i . One arrives, thus, at a ther- 
modynamic Landau-type description of the first-order phase transitions. Notice, 
that the effective thermodynamic description appears due to the assumption that 
H (p, q) = and, therefore, that all the processes take an infinitely long time. 



5. Fermions 

5.1. Free fermion Keldysh action 

Consider a single quantum state, with the energy e - This state is populated 
by spin-less fermions (particles obeying the Pauli exclusion principle). In fact, 
one may have either zero, or one particle in this state. The secondary quantized 
Hamiltonian of such a system has the form: 

H=e Jc, (5.1) 

where c< and c are fermion creation and annihilation operators of the state eo. 
They obey standard anf/commutation relations: {c,a}+ = 1 and {c,c} + = 
{c^ , c'}+ — 0, where { , }+ stands for the anti-commutator. 

One can now consider the evolution operator along the Keldysh contour, C 
and the corresponding "partition function", Z = 1, defined in exactly the same 
manner as for bosonic systems: Eq. i2.2\ . The trace of the equilibrium density 

5 Indeed, Tp = d t q = -DV^q + qU£' q = T(-DV%p + pU£' q ). This non-trivial fact reflects 

the existence of an accidental conservation law: H(p(t; r), q(t; r)J = - locally] While from the 
general principles only the total global energy has to be conserved. 
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matrix is Tr{po} = 1 + p(^o), where the two terms stand for the empty and the 
singly occupied state. One divides the Keldysh contour into (27V — 2) time in- 
tervals of length St ~ 1/JV — ► and introduces resolutions of unity in 2N points 
along C, Fig. II. 21 . The only difference from the bosonic case in section |2~TI is 
that now one uses a resolution of unity in the fermionic coherent state basis [17]: 



(5.2) 



where ipj and ipj are mutually independent Grassmann variables. The rest of the 
algebra goes through exactly as in the bosonic case, section l2~T1 As a result, one 
arrives at: 



z = 5vM/yn[**]e"'" 



(5.3) 



where the 2N x 27V matrix Q -} stands for: 



^33 



-1 


-p( £ o) " 


1-h -1 




1-h -1 




1 


-1 




1+h -1 




1+h -1 . 



(5.4) 



and h = ieoSt- The only difference from the bosonic case is the negative sign 
before the p(eo) matrix element, originating from the minus sign in the (— tp2N\ 



Reminder: the fermionic coherent state \ip) = (1 — ipc')\0), parameterized by a 
Grassmann number ip (such that {ijj, = {ip, c} + = 0), is an eigenstate of 
the annihilation operator: c\tp) = rp\ip). Similarly: (ip\c> — where ^ is 

another Grassmann number, unrelated to tp. The matrix elements of a normally 
ordered operator, such as e.g. the Hamiltonian, take the form 
(ip\H (c' , c)\ip') = H(ip,ip')(i(>\i(>'). The overlap between any two coherent 
states is (ip\ip') = 1 + ijjip' = exp{^V'}- The trace of an operator, A, is 
calculated as: Tr{j4} = JJ dip dip e~^(— ip\A\ip), where the Grassmann 
integrals are defined as: J dip 1 = and Jdipip = 1. 



42 



A. Kamenev 



coherent state in the expression for the fermionic trace. To check the normaliza- 
tion, let us evaluate the determinant of such a matrix: 

det^- 1 ] = l + p( £o )(l-/ l 2 ) Ar - 1 « l+pMe^) 2 ^" 1 ) -> l+p(e ).(5.5) 

Employing the fact that the fermionic Gaussian integral is given by the deter- 
minant (unlike the inverse determinant for bosons) of the correlation matrix, 
|Appendix A| one finds: 

det^er 1 ! 

Z = I , J = 1 , (5.6) 
Tr{p } 

as it should be. Once again, the upper-right element of the discrete matrix, 
Eq. ( 15. 4t . is crucial to maintain the correct normalization. 

Taking the limit N — > oo and introducing the continuous notations, tpj — > 
tp(t), one obtains: 

Z = jv^ij) e iS &<^ = Jv^ exp ji ^ [^(^"VM] dt J , (5.7) 
where according to Eqs. d5.3l > and ( 15. 4> the action is given by 

2N 



3=2 



*^3 fc. eo^iVj-i 



+ i -01 ("01 + p(eo)V>2Jv) , (5.8) 



where <Jtj = tj — tj-i = ±<5{. Thus the continuous form of the operator 1 is 
the same as for bosons, Eq. (12. 1 31 : Q^ 1 = idt — £o- Again the upper-right ele- 
ment of the discrete matrix (the last term in Eq. ( 15. 8t ). that contains information 
about the distribution function, is seemingly absent in the continuous notations. 

Splitting the Grassmann field ip(t) into the two components tp+{t) and ip-(t) 
that reside on the forward and the backward parts of the time contour correspond- 
ingly, one may rewrite the action as: 



S= dtil>+{t)[id t -e Q ]il> + {t)- dtil>-(t)[idt-eoty-{t), (5.9) 



where the dynamics of ip + and ip_ are not independent from each other, due to 
the presence of non-zero off-diagonal blocks in the discrete matrix, Eq. \5A\ . 
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The four fermionic Greens functions: CJ T ( T ) and Q <ty> ^ are defined in the 
same way as their bosonic counterparts, Eq. \2.\6\ : 

WW- O = J y M = rr f - i = -n F e ru > ; 

Tr{po} 

a>-w + (t')) = iQ > m = Tr{c ^ } ,)/3o} - (1-^)^-^(5.10) 

(V+W^+(0> = iG T (t,t')=6(t-t')ig > (t,t')+6(t' -t)ig<(t,t'); 
(^_(t)$_(0) = iG f {tJ)=e{1?-t)iQ > {tJ) + 9{t-1?)iQ < {tJ); 

The difference is in the minus sign in the expression for due to the anti- 
commutation relations, and the Planck occupation number is exchanged for the 
Fermi one: n — > rip = p(e )/(l + p(eo)). Equations ( I2. 17t and ( I2. 1 81 hold for 
the fermionic Green functions as well. 

5.2. Keldysh rotation 

It is customary to perform the Keldysh rotation in the fermionic case in a different 
manner from the bosonic one. Define the new fields as: 

1 p 1 (t) = ^=(i>+(t)+^-(t)); 1 p 2 (t) = ^=(i>+(t)-^-(t)). (5.11) 

This line is exactly parallel to the bosonic one, Eq. J2. 19b . However, following 
Larkin and Ovchinnikov [20], it is agreed that the "bar" fields transform in a 
different way: 

&(f) = -i= ($+(<)- $_(«)); &(t) = ^$+(t)+^-(*))- (5.12) 

The point is that the Grassmann fields if) are not conjugated to if), but rather are 
completely independent fields, that may be chosen to transform in an arbitrary 
manner (as long as the transformation matrix has a non-zero determinant). No- 
tice, that there is no issue regarding the convergence of the integrals, since the 
Grassmann integrals are always convergent. We also avoid the subscripts cl and 
q, because the Grassmann variables never have a classical meaning. Indeed, one 
can never write a saddle-point or any other equation in terms of tfj, ip rather they 
must always be integrated out in some stage of the calculations. 
Employing Eqs. ( 15.1 It . ( 15.12b along with Eq. ( 15.101 1. one finds: 

- i{Mt)Mt')) ^G ab =( QR{t Q gl^] ) » (5-13) 
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where hereafter a,b = (1, 2). The presence of zero in the (2, 1) element of this 
matrix is a manifestation of identity (12. 1 8i . The retarded, advanced and Keldysh 
components of the Green function are expressed in terms of ^ T ( T ) and in 
exactly the same way as their bosonic analogs, Eq. J2.2H . and therefore posses 
the same symmetry properties: Eqs. J2.22t - ll2.24t , An important consequence of 
Eqs. d2~23l . d2~24l is: 

Tr{a 1 oS 2 o...o4}(t,i) = 0, (5.14) 

where the circular multiplication sign involves integration over the intermediate 
times along with the 2x2 matrix multiplication. The argument (t, t) states that 
the first time argument of Q\ and the last argument of Qu are the same. 

Notice that the fermionic Green function has a different structure from its 
bosonic counterpart, Eq. 12.201 : the positions of the R, A and K components in 
the matrix are exchanged. The reason, of course, is the different convention for 
transformation of the "bar" fields. One could choose the fermionic convention 
to be the same as the bosonic (but not the other way around!), thus having the 
same structure, Eq. J2.20i . for fermions as for bosons. The rationale for the 
Larkin-Ovchinnikov choice, Eq. J5.13> . is that the inverse Green function, C? _1 
and fermionic self energy £ p have the same appearance as Q: 

whereas in the case of bosons G~ x , Eq. ( 12.291 1. and S, Eq. (13. lOi . look differently 
from G, Eq. (I2.20> . This fact gives the form Eq. J5.13> . d5.15> a certain technical 
advantage. 

For the single fermionic state (see. Eq. J5.10» : 

-i6{t - t') £>-«<>(*-*') -> ( e - e + iQ)- 1 ; 
i6{t' - t) e~ iea ( *~*' ) -► (e - e - iO^ 1 ; (5.16) 
- 2n F ) e' 1 ^-^ (1 - 2njr(e))(-27r»)<y(e - e ) . 

where the r.h.s. provides also the Fourier transforms. In thermal equilibrium, one 
obtains: 

^(e)=tanh^(^(e)-^(e)) . (5.17) 

This is the FDT for fermions. As in the case of bosons, the FDT statement 
is a generic feature of an equilibrium system, not restricted to the toy model. 



G R (t,t>) = 
Q A (t,t') = 
G K (t,t>) = 
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In general, it is convenient to parameterize the anti-Hermitian Keldysh Green 
function by a Hermitian matrix T = as: 

Q K =Q R oT-ToQ A , (5.18) 

The Wigner transform of T{t, t') plays the role of the fermionic distribution func- 
tion. 

One may continue now to a system with many degrees of freedom, counted by 
an index k. To this end, one simply changes: eo — > efe and perform summations 
over k. If A: is a momentum and e/- = k 2 / (2m), it is instructive to transform to the 
real space representation: i/j(t;k) — > ip(t;r) and = k 2 /{2m) = — (2m)~ 1 V 2 . 
Finally, the Keldysh action for a non-interacting gas of fermions takes the form: 

S $,il>] = ffdxdx'J2 ^ a {x)[g- l {x,x')) ab ^ b {x') , (5.19) 

a.b=l 

where x — (i; r) and the matrix correlator [£ _1 ] a i> has the structure of Eq. 15. 151 
with 



[gXW&x')]- 1 =5(x-x') (id t , + ^-\7 2 r ) 

\ 2m J 



(5.20) 



Although in continuous notations the R and the A components look seemingly 
the same, one has to remember that in the discrete time representation, they are 
matrices with the structure below and above the main diagonal correspondingly. 
The Keldysh component is a pure regularization, in the sense that it does not 
have a continuum limit (the self-energy Keldysh component does have a non- 
zero continuum representation). All this information is already properly taken 
into account, however, in the structure of the Green function, Eq. 15.131 . 

5.3. External fields and sources 

Let us introduce an external time-dependent scalar potential —V(t) defined along 
the contour. It interacts with the fermions as: Sy = f c dtV{t)^>{t)if>(t). Ex- 
pressing it via the field components, one finds: 



S v =ldt[V+$+il>+-V-$-il>J\=jdt[V cl ($ + 1> + ^ 

UtlVai^x+^^+VSi^+^i)] , (5.21) 



-oo 

oo 
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where the V c i and the V q components are defined in the standard way for real 
bosonic fields: V c u q \ = (V+±V_)/2. Notice that the physical fermionic density 
(symmetrized over the two branches of the contour): g = | i'4 > +' l P++'4 > -' l P—) is 
coupled to the quantum component of the source field, V q . On the other hand, 
the classical source component, V c i, is nothing but an external physical scalar 
potential, the same at the two branches. 

Notations may be substantially compactified by introducing vertex gamma- 
matrices: 

**= ( 1 V *fl= 1 " ' 



1=[ Q J ) ; ^={1 o ] (5-22) 

With the help of these definitions, the source action 15.211 may be written as: 

oo 2 00 

S V = Jdt^2 [Vci $alab<Pb + V q i>allM = U ^„f^ , (5-23) 
-oo Q ' b=1 -oo 

where the summation index is a = (cl, q). 

Let us define now the "generating" function as: 

Z[V cll V q ] = (e iSv ) , (5.24) 

where the angular brackets denote the functional integration over the Grassmann 
fields 'ip and ip with the weight of e lS ° specified by the fermionic action d5 . 1 91 . 
In the absence of the quantum component, V q = 0, the source field is the same 
at both branches of the contour. Therefore, the evolution along the contour still 
brings the system back to its exact initial state. Thus one expects that the classical 
component alone does not change the fundamental normalization, Z = 1 . As a 
result: 

Z[V cl ,V g = 0] = l. (5.25) 

One may verify this statement explicitly by expanding the action in powers of 
Vci and employing the Wick theorem. For example, in the first order one finds: 
Z[Vd, 0] = 1 + / dt Vd(t) Tr{Y l G(t, t)} = 1, where one uses that j cl = 1 along 
with Eq. J5. 14b . It is straightforward to see that for exactly the same reasons all 
higher order terms in V c i vanish as well. 

A lesson from Eq. J5.25t is that one necessarily has to introduce quantum 
sources (that change sign between the forward and the backward branches of the 
contour). The presence of such source fields explicitly violates causality, and thus 
changes the generating function. On the other hand, these fields usually do not 
have a physical meaning and play only an auxiliary role. In most cases one uses 
them only to generate observables by an appropriate differentiation. Indeed, as 
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was mentioned above, the physical density is coupled to the quantum component 
of the source. In the end, one takes the quantum sources to be zero, restoring the 
causality of the action. Notice that the classical component, V c u does not have to 
be taken to zero. 

Let us see how it works. Suppose we are interested in the average fermion 
density g at time t in the presence of a certain physical scalar potential V c i(t). 
According to Eqs. (I5.2H and 15 . 241 it is given by: 



g(t; V cl 



2 SV q (t) 



Z[V ch V q ] 



v„=o 



(5.26) 



The problem is simplified if the external field, V c i, is weak in some sense. Then 
one may restrict oneself to the linear response, by defining the susceptibility: 



n R (t,t') 



W^) e[t]Vcl) 



V r ,=0 



1 S 2 Z[V ch V q ] 

2 SV d (t')SV g {t) 



.(5.27) 



v a =v cl =o 



We add the subscript R anticipating on physical grounds that the response func- 
tion must be retarded (causality). We shall demonstrate it momentarily. First, let 
us introduce the polarization matrix as: 



IT a/3 (M') 



5 2 lnZ[V] 



2 SV (t')SV a {t) 





n fl (M') 



n A (t,t') 
u K (t,f) 



(5.28) 



Due to the fundamental normalization, Eq. ( I5.25t . the logarithm is redundant for 
the R and the A components and therefore the two definitions J5.27t and J5.28i 
are not in contradiction. The fact that H cLcl = is obvious from Eq. J5.25> . 
To evaluate the polarization matrix, II, consider the Gaussian action, Eq. ( 15. 191 . 
Adding the source term. Eq. J5.23> . one finds: Sq + Sv = J dt 'tpiG -1 + V a A f a ]i[> '. 
Integrating out the fermion fields %f>, ip according to the rules of fermionic Gaus- 
sian integration, |Appendix A| one obtains: 

Z[V]^^dct{tg- 1 +iV a r}=dct{l + gV a r}=e Tnn{1+dVar \ (5.29) 

where one used normalization, Eq. d5.6i . Notice, that the normalization is exactly 
right, sinceZ[0] = 1. One may now expand ln(l+£? V Q 7°) to the second order in 
V and then differentiate twice. As a result, one finds for the polarization matrix: 



n 



{t,t') = - l -Tv{rG{t,t')^Q{t',t)) 



(5.30) 



48 



A. Kamenev 



Substituting the explicit form of the gamma-matrices, Eq. d5.22> . and the Green 
functions, Eq. (I5.13> . one obtains for the response and the correlation compo- 
nents: 



gR(A) ^ t ')gK( t i 7t) + QK (tj tl) gA{R) y ^ f) 



(5.31) 



n K (t, t')=- 1 - [g K (t, t')g K (t', t)+g R (t, t')g A (t>, t)+g A (t, t')g R (t>, t)] . 

From the first line it is obvious that H R W (t, t') is indeed a lower (upper) trian- 
gular matrix in the time domain, justifying their superscripts. Moreover, from the 
symmetry properties of the fermionic Green functions (same as Eq. (I2.22» one 
finds: H R = [II y and II A = — [II A ]^. As a result, the polarization matrix, ft, 
possesses all the symmetry properties of the bosonic self-energy fi, Eq. (13. 101 . 



Exercise: In the stationary case: Q(t, t') = Q(t — t'). Fourier transform to the energy domain and 
write down expressions for Ti(ui). Assume thermal equilibrium and, using Eq. 15. 171 . rewrite your 
results in terms of Q R ( A ^ and the equilibrium distribution function. Show that in equilibrium, the 
response, H R ( A ) (u>), and the correlation, H K (lu), functions are related by the bosonic FDT: 



H = coth^(n B H-ir>)) . 



(5.32) 



Equation J5 .3 1 1 for H R constitutes the Kubo formula for the density-density 
response function. In equilibrium it may be derived using the Matsubara tech- 
nique. The Matsubara routine involves, however, the analytical continuation from 
discrete imaginary frequency uj m to real frequency ui. This procedure may prove 
to be cumbersome in specific applications. The purpose of the above discussion 
is to demonstrate how the linear response problems may be compactly formu- 
lated in the Keldysh language. The latter allows to circumvent the analytical 
continuation and yields results directly in the real frequency domain. 

5.4. Tunnelling current 

As a simple application of the technique, let us derive the expression for the 
tunnelling conductance. Our starting point is the tunnelling Hamiltonian: 



H = E [4 C) 4c* + A d) Adk] + £ \Tw Ad k , + T* kk , d\,c k 



(5.33) 



where the operators c k and d k > describe fermions in the left and right leads, while 
T kk > are tunnelling matrix elements between the two. The current operator is: 

j =TtY<k C W = Efe c\ c k]- = -i Efcfc' T kk> c{d k > - T£ fc , d\,c k . 



Many-body theory of non-equilibrium systems 



49 



To describe the system in the Keldysh formalism, one introduces the the four- 
component spinor: ip k = (M>f k \ "tp^k > $ik ' $2k)> a similarly one f° r the fields 
without the bar, and the 4x4 matrices: 

Qk= [ o jjwj. o y Jkk, -[-iT^ o j- (5 - 34) 

In addition to the already familiar Keldysh structure the spinors and matrices 
above possess the structure of the left-right space. In terms of these objects the 
action and the current operator take the form: 

oo 

S=jdt'S^\l> k Skk'Gk 1 - Tk,w 4>k> ; J(t) = - V] V^Jfc.fc/Vv ■ (5.35) 



The current is expressed through the 7 9 vertex matrix in the Keldysh space be- 
cause any observable is generated by differentiation over the quantum compo- 
nent of the source field (the classical component of the source does not change 
the normalization, Eq. (I5.25» . 

One is now in a position to calculate the average tunnelling current up to the 
second order in the matrix elements Tk,k' ■ To this end one expands the action up 
to the first order in T^w, and applies the Wick theorem: 

oo 

J{t)=iJdt l Y J ^{hkiGk-{t,t l )T k , k g k {t\t)} (5.36) 

— oo kk> 

oo 

= jdt'Y, \Tkk> i 2 Tr {rg { k c) a, f)r l g { f (f, t) - *p&*> a, t>)r l g { k c) v, *) } = 

-oo kk ' 
oo 

/ dt'y^jTkk 

J i, i, I 



kk i 



r (c)R r (d)K r (c)K r (d)A _ r (d)R r (c)K _ r (d)K r (c)A 
y k(t,t') y k'(t',t) +y k{t,t') y k'(t',t) y k'{t,t') y k{t',t) y k'(t,t') y k(t',t) 



Now let us assume a stationary situation so that all the Green functions depend on 
the time difference only. We shall also assume that each lead is in local thermal 
equilibrium and thus its Green functions are related to each other via the FDT: 

G k c)K (e) = {l-2n ( p ) (e))[g k c)R (e)-g i k c)A {e)}. Similarly for the "d"-lead with 

a different occupation function rip\e). As a result, one finds for the tunnelling 
current: 

oo 

J=\- [»£>(«) - 4\e)]J2\T kk , | 2 - Q\? A ] [Q^ - G^ A ) . (5.37) 

_"L n kk' 
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If the current matrix elements may be considered as approximately momentum 
independent: |Tfefe'| 2 |T| 2 , the last expression is reduced to: 



oo 

J = Att\T\ 2 fde [nP(e) - nf{e)} i»(e) i/W(e) , (5.38) 



where the density of states (DOS) is defined as: 



5.5. Interactions 

Consider a liquid of fermions that interact through instantaneous density-density 
interactions: Hint = ~^ffdrdr' : g(r)U(r — r')g{r') :, where g(r) = c\c r 
is the local density operator and : . . . : stands for normal ordering. The corre- 
sponding Keldysh contour action has the form: Si n t = \ JcdtJJdrdr'U(r — 
r') , 4) r '4i r iip r il) r i . One may now perform the Hubbard-Stratonovich transforma- 
tion with the help of a real boson field ip(t; r), defined along the contour: 

^ J dtJJ drdr' U(r—r')rp r tp r r i^ r tp r / p ij dt^JJ drdr' ip r U~) (p r i 4- Jdrip r tp r rp r ^ 

e ° = Vipe c ,(5.40) 

where U^ 1 is a kernel, that is inverse to the interaction potential: U^ 1 o U = 
1. One notices that the auxiliary bosonic field, ip, enters the fermionic action 
in exactly the same manner as a scalar source field. Following Eq. ( I5.21l i. one 
introduces (p c i( q ) = (tp+ ± ¥?— )/2 and rewrites the fermion-boson interaction 
term as 'i'afalab^b > wnere summations are assumed over a,b = (1,2) and a = 
(cl, q). The free bosonic term takes the form of: ^ipU^ 1 ^ — > (pallia™ <pp. 

At this stage the fermionic action is Gaussian and one may integrate out the 
Grassmann variables in the same way it was done in Eq. ( 15.291 1. As a result, one 
finds for the generating function, Eq. J5.24I . of the interacting fermionic liquid: 

„ i J dt JJdrdr'ip U^aiip + Tr In [l+g (V a +ip a )y a ] 

Z[V] = \Vip e —> . (5.41) 

Quite generally, thus, one may reduce an interacting fermionic problem to a the- 
ory of an effective non-linear bosonic field (longitudinal photons). Let us demon- 
strate that this bosonic theory possesses the causality structure. To this end, one 
formally expands the logarithm on the r.h.s. of Eq. ( 15. 4H . Employing Eq. ( 15 . 1 4i 
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and recalling that A f cl = 1, one notices that for (p q = V q = the bosonic action 
is zero. As a result, Eq. J2.32t holds. 

To proceed we shall restrict ourselves to the, so called, random phase approx- 
imation (RPA). It neglects all terms in the expansion of the logarithm beyond 
the second order. The second order term in the expansion is conveniently ex- 
pressed through the (bare) polarization matrix IP <" (see Eq. J5.30» of the non- 
interacting fermions. The resulting effective bosonic theory is Gaussian with the 
action: 



S R PA[<f,V]=ljdtdt'JJdrdr' <p (u^ai-Uj <p -2<pUy-VUV 



(5.42) 



One notices that the bare polarization matrix plays exactly the same role as of the 
self-energy, S, cf. Eqs. J3.9t >. ( 13.1 Oi l, in the effective bosonic theory. As a result, 
the full bosonic correlator ([/ _1 (7i — II) possesses all the causality properties, 
listed in section l2~4l 

Finally, let us evaluate the dressed polarization matrix of the interacting fermi- 
liquid in the RPA. To this end one may perform the bosonic Gaussian integra- 
tion in the RPA action J5.42t to find the logarithm of the generating function: 

i In Z RPA [V} = v(tt + fl(U- 1 a 1 - ft)" 1 !!) V. Finally, employing the defi- 
nition of the polarization matrix, Eq. J5.28l >. and performing simple matrix alge- 
bra, one finds: 

ti-RPA =fto (l-^Uofiy 1 . (5.43) 

It is straightforward to demonstrate that the dressed polarization matrix possesses 
the same causality structure as the bare one, Eq. J5.28b . For the response com- 
ponent of the dressed polarization, Tl R , pA , the second factor on the r.h.s. of 
Eq. J5.43l > may be considered as a modification of the applied field, V c i ■ Indeed, 
cf. Eq. J5.27> . q = Tl^ pA o V c i — Yl R o V£ cr , where the screened external po- 
tential V c f r is given by: V c f r = (l - a x U o fl R ^j oV ci . This is the RPA result 
for the screening of an external scalar potential. 



5.6. Kinetic equation 

According to Eq. J5.3QI > to evaluate the bare (and thus RPA dressed, Eq. ( I5.43t ) 
polarization matrix, one needs to know the fermionic Green function, Q. While 
it is known in equilibrium, it has to be determined self-consistently in an out- 
of-equilibrium situation. To this end one employs the same idea that was used 
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in the bosonic theory of chapter[3] Namely, one writes down the Dyson equation 
for the dressed fermionic Green function: 

(g^ 1 -E f ) o£ = i, (5.44) 

where the subscript "0" indicates the bare Green function. The fermionic self- 
energy, Si? turns out to have the same structure as Q~ x , Eq. (15. 151 . Thus the R 
and A components of the Dyson equation take a simple form: 

^ + _L Vr 2 ) = e£ (A) o e R(A) . (5.45) 

Employing the parameterization Q K = Q Tl oT—J-oQ A , where J 7 is aHermitian 
matrix, along with Eq. j5.45L one rewrites the Keldysh component of the Dyson 
equation as: 

= Ef - (Sp o T - T o E£) - -i I col [F] . (5.46) 

This equation is the quantum kinetic equation for the distribution matrix T. Its 
l.h.s. is the kinetic term, while the r.h.s. is the collision integral with E^ having 
the meaning of an "incoming" term and E^ o T ' — T ' o E^ that of an "outgoing" 
term. 

The simplest diagram for the fermionic self-energy matrix, E^ fc , is obtained 
by expanding the Hubbard-Stratonovich transformed action, Eq. J5.40i . to the 

second order in the fermion-boson interaction vertex, $ a <Pa A /%b'tpb> ar, d applying 
the Wick theorem for both fermion and boson fields. As a result, one finds: 

Sf(M') = (laJ cd (t,t')f db ) ( Va (t)Mt'))=(r l G(t,t>)rfiD K (t,ti); 

+ (f G(t, if) j q J b iD R (t, t')+(y q S(t, t') rfiD A (t, 0, (5.47) 

where summations over all repeated indexes are understood and the spatial argu- 
ments have the same general structure as the time ones. The boson Green func- 
tion is denoted as ((p a (t)<pp(t')) = iD al3 (t, if). Finally one finds for the R, A 
(i.e. (1, 1) and (2, 2)) and K (i.e (1,2)) components of the fermionic self-energy: 

T,f A \%t')=i{g R ^\t,t')D K {t,t') + g K {t,t')D R ^ (t,f)) ; (5.48) 

s£(t, if) = i(g K (t, t')D K {t, t') + g R (t, t')D R {t, t') + g A (t, t')D A (t, t')) 
- i(g K (t, t')D K (t', t) + (g R (t, t') - g A (t, t')) (D R (t, t') - D A (t, t'))) , 
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where in the last equality one had used that Q R W (t, t')D A{ * R:) (t, t') = 0, since 
these expressions have no support in the time domain (see, however, the foot- 
note in section [HJ For the same reason: T,f(t,t') = i(g A (t,t')D R (t,t') + 
Q R (t,t')D A (t,t')^j = 0. As expected, the retarded and advanced components are 
lower and upper triangular matrices correspondingly, with T, R — [S A ]\ while 
Yi K = — [S 7 ^. Notice the close resemblance of expressions J5.48> to their 
bosonic counterparts, Eqs. (13.1 4i — (13 . 1 6i . 

If one understands the bosonic Green function, D, as the bare instantaneous 
interaction potential (i.e. D R = D A = U(r - r')S(t - t') and D K = 0), one 
finds: = = iUg K {t, t)S(t - t') and Ef = 0. In this approximation the 
r.h.s. of the kinetic equation (I5.46> vanishes (since T is a symmetric matrix) and 
so there is no collisional relaxation. Thus one has to employ an approximation 
for D that contains some retardation. The simplest and most convenient one is 
the RPA, where D = {U^&i - ft) -1 , cf. Eq. (I5.42i . with a matrix ft that is 
non-local in time. This relation may be rewritten as the Dyson equation for D, 
namely (U^ 1 <ti — II) o D = 1. One may easily solve it for the three components 
of D and write them in the following way: 



Performing the Wigner transform following sections 1331 13. 61 the kinetic term 
(the l.h.s. of Eq. J5.46» is exactly the same as for the complex boson case (one 
has to take into account the gradient terms to obtain a non-zero result for the WT 
of the commutator). The result is (cf. Eq. J3.17» : 



where Vk = d^k, E is an external electric field and the collision integral, I co i, 
is i times the WT of the r.h.s. of Eq. (15.461 . On the r.h.s. one may keep only the 
leading terms (without the gradients). One also employs a parameterization of the 
Keldysh component of the fermionic Green function through the corresponding 
distribution function: Q K — > g K — fF(g fl — g A ), where fp(T, p, k) is the WT 
of T. Assuming, for brevity, a spatially uniform and momentum isotropic case, 
one may restrict oneself to fp(r, e&) = fp(r, e). As a result, one finds for the 
collision integral: 






(5.50) 




i 



(5.51) 



x[(n 



IT 4 ) (1 - f F (e - w)fr(e)) - n K (f F (e) - f F ( 



e-w))], 
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where n a/3 = H a ^(uj, q), while the time index, r, is suppressed for brevity and 
the notation 

A g (e, k) = ±- (g R (e, k) - g A (e, k)) (5.52) 

is introduced. For free fermions A g (e, k) = S(e — ef.). At this stage one may ob- 
serve that if the bosonic system is at equilibrium: II A = coth(w /2T) [IT^ — II 1 , 
then the fermionic collision integral is nullified by: 

f F (e) = tanh ^ . (5.53) 

Indeed, 1 — tanh(fe— a) tanh(&) = coth(o)(tanh(6)— tanh(£>— a)). One should 
take into account, however, that the bosonic degrees of freedom are not indepen- 
dent from the fermionic ones. Namely, components of the polarization matrix II 
are expressed through the fermionic Green functions according to Eq. J5.3 II . In 
the WT representation these relations take the form: 

T1 R - IT 4 = m /de'5^A,(e', k')A e (e' -u,k'- q) [f F (e' - w) - f F (e')] ; 
^ w 

H K (uj 7 q) = i 7 r/'de'^Ag(e' ! fc , )A g (e'-w,fc'-g) [f F (e'-w)f F (e') - 1] .(5.54) 

Due to the same trigonometric identity the equilibrium argument can be made 
self-consistent: if the fermionic system is in equilibrium, Eq. ( I5.53t . then com- 
ponents of II satisfy the bosonic FDT, Eq. P.22> . 

One may substitute now Eqs. J5.54i into Eq. ( I5.5H to write down the closed 
kinetic equation for the fermionic distribution function. Most conveniently it is 
done in terms of the occupation numbers, defined as f F = 1 — 2n 6 : 



^=H<U,I,'\I 



rv ( l-n e ) ( 1 - ) - n e n e ^ ( l-n e / ) ( 1-n^) 
where the transition probability is given by: 



,(5.55) 



M(e, lu) =4vr^ \D r {uj, q)\ 2 A g (e-w, k-q)A g (e', fc')A g (e'-w, k'-q). (5. 56) 

Equation d5.55l > is a generic kinetic equation with a "four-fermion" collisional 
relaxation. The first term in the square brackets on its r.h.s. may be identified as 
"in", while the second one as "out". Each of these terms consists of the prod- 
uct of four occupation numbers, giving a probability of having two initial states 

6 To derive this expression one should add and subtract n e n e /_ [J n (! / n e —^ in the square brackets. 
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Fig. 7. Structure of the four-fermion collision integral. The full lines are fermionic Green func- 
tions; the wavy lines are the RPA screened interaction potential. The fermionic loop represents the 
polarization matrix, II (id, <?). 



occupied and two final states empty. For n(e) given by the Fermi function the 
"in" and the "out" terms cancel each other. Therefore in thermal equilibrium the 
components of the dressed fermionic Green function must satisfy the FDT: 

£ K = tanh^(g R -C? A ). (5.57) 

The structure of the transmission probability M is illustrated in Fig. 15.61 The 
three factors of A g enforce that all three intermediate fermionic particles must 
satisfy energy-momentum conservation (stand on mass-shell), up to the quasi- 
particle life-time. The real factor |D fl | 2 is associated with the square of the 
matrix element of the screened interaction potential (in the RPA). 



6. Disordered fermionic systems 

6.1. Disorder averaging 

We consider fermions in the field of a static (quenched) space-dependent scalar 
potential Udi S (r). The potential is meant to model the effect of random static 
impurities, dislocations, etc. Since one does not know the exact form of the 
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potential, the best one can hope for is to evaluate the statistical properties of var- 
ious observables, assuming some statistics for Udis{r). It is usually a reasonable 
guess to prescribe a Gaussian distribution for the potential. Namely, one assumes 
that the relative probability for a realization of the potential to appear in nature is 
given by: 



where v is the bare fermionic DOS at the Fermi level and r, called the mean-free 
time, measures the strength of the random potential. 

In this chapter we concentrate on non-interacting fermions. We would like 
to evaluate, say, the response function, H R , in presence of the random potential 
and average it over the realizations of Udis with the weight given by Eq. ( 16. 1> . 
The crucial observation is that the response function, 11^, may be defined as 
variation of the generating function, Eq. J5.27l >. and not the logarithm of the 
generating function. More precisely, the two definitions with, Eq. J5.28> . and 
without, Eq. J5.27t . the logarithm coincide due to the fundamental normalization, 
Eq. J5.25b . This is not the case in the equilibrium formalism, where the presence 
of the logarithm (leading to the factor Z^ 1 after differentiation) is unavoidable 
in order to have the correct normalization. Such a factor Z^ 1 = Z~ 1 [Udi S ] 
formidably complicates the averaging over Udis- Two techniques were invented 
to perform the averaging: the replica trick [21] and the super-symmetry (SUSY) 
[22]. The first one utilizes the observation that InZ = \im n ^o(Z n — to 
perform calculations for an integer number, n, replicas of the same system and 
take n — > at the end of the calculations. The second one is based on the fact that 
Z^ 1 of the non-interacting fermionic system equals to Z of a bosonic system in 
the same random potential. One thus introduces an additional bosonic replica 
of the fermionic system at hand. Both of these ideas have serious drawbacks: 
the replica technique requires analytical continuation, while the SUSY is not 
applicable to interacting systems. 

The Keldysh formalism provides an alternative to these two methods by insur- 
ing that Z = 1 by construction. One may thus directly perform the averaging of 
the generating function, Eq. ( I5.24t . over realizations of Udis- Since the disorder 
potential possesses only the classical component (it is exactly the same on both 
branches of the contour), it is coupled only to j cl = 1. The disorder-dependent 
term in the averaged generating function has the form: 



V[U dis ] 



e 



-JdrU 2 dts (r) 



(6.1) 




Jdr TTVTU 2 dis {r)-iU dts ,{r) J dt$ t 7 cl $ t 



(6.2) 
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-4^r f dr JJ dtdt ' (^ a V>?)«/V> t 6 ,) 
= e , 

where a, 6 = 1,2, and there is a summation over repeated indexes. One can 
rearrange the expression in the exponent on the r.h.s. of the last equation as 
(fyipt) (V't'V't') = - (i't'^t) 7 and then use the Hubbard-Stratonovich 

matrix field, Q = Qf t ,(r): 



-fdr 



= VQ 



^Tr{QoQ}-^ JJ dtdt'Qf^\,r t 



J, (6.3) 



where the spatial coordinate, 7-, is suppressed in both Q and ip. At this stage the 
average action becomes quadratic in the Grassmann variables and they may be 
integrated out leading to the determinant of the corresponding quadratic form: 
Q^ 1 + V a j a + j^Q- All the matrices here should be understood as having a 
2x2 Keldysh structure along with an N x N structure in discrete time. One thus 
finds for the disorder averaged generating function: 



Z[V]= VQe 



iS[Q;V] 



where 



iS[Q-V] 



At 



drTr{Q 2 } + Trm 



i 



(6.4) 



(6.5) 



As a result, one has traded the initial functional integral over the static field 
Udis(r) for the functional integral over the dynamic matrix field Qt.f(r). At 
first glance, it does not strike as a terribly bright idea. Nevertheless, there is a 
great simplification hidden in this procedure. The point is that the disorder po- 
tential, being ^-correlated, is a rapidly oscillating function. On the other hand, as 
one will see below, the Q-matrix field is a slow (both in space and time) function. 
Thus it represents the true macroscopic (or hydrodynamic) degrees of freedom 
of the system, that happen to be the diffusively propagating modes. 



6.2. Non-linear a-model 

To execute this program, one first looks for the stationary configurations of the 
action d6.5i . Taking the variation over Qt,# (r), one obtains: 



,(r) 



i 



2t 



Q 



(6.6) 



7 The minus sign originates from commuting the Grassmann numbers. 
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where Q denotes a stationary configuration of the fluctuating field Q. For the 
purpose of finding the stationary configurations one has omitted the small source 
field, V. It is important to notice that the spatially non-local operator [Gq + 

j-Q] 1 (t, t'; r, r') on the r.h.s. is taken at coinciding spatial points r' = r. 

The strategy is to find first a spatially uniform and time-translationally invari- 
ant solution of Eq. J6.6t : Q t t ,, and then consider space and time-dependent 
deviations from such a solution. This strategy is adopted from the theory of 
magnetic systems, where one first finds the uniform static magnetized configu- 
rations and then treats spin-waves as smooth perturbations on top of such static 
uniform solutions. From the structure of Eq. ( 16. 6> one expects that the saddle- 
point configuration Q possesses the same structure as the fermionic self-energy, 
Eq. ( 15 . 1 5I > (more accurately, one expects that among the possible saddle points 
there is a "classical" one, that satisfies the causality structure, Eq. d5.15» . One 
looks, therefore, for a solution of Eq. ( 16. 6> in the form of: 



Q c = A £ = 



Af 




Af 
A* 



(6.7) 



Substituting this expression in Eq. ( 16. 6t . one finds 



a5^=— ■ 



R(A) 



--S-A 

2t 



R(A) 



2t 



A 



R{A) 



= ±1, (6.8) 



where ^ = k 2 / (2m) — /i and one adopts J^k • ■ ■ = v I d£k ■ ■ •> where v is the 
DOS at the Fermi surface. The summation over momentum appears because the 
matrix on the r.h.s. is taken at coinciding spatial points. The signs are chosen so 
as to respect causality: the retarded (advanced) Green function is analytic in the 
entire upper (lower) half-plane of complex energy e. One has also assumed that 
1 / (2t) <C /i. The Keldysh component, as always, may be parameterized through 
a Hermitian distribution function matrix: A K = k R oT — Tok A = 2T t , where 
the distribution function T t is not fixed by the saddle point equation ( 16. 6i and 
must be determined through the boundary conditions. As a result one obtains: 



A e = 



1 2T t 
-1 



(6.9) 



Transforming back to the time representation, one obtains A 



R(A) _ 



±5(t -t'T 



0), where =f0 indicates that the (5-function is shifted below (above) the main 
diagonal, t = t' . As a result, Tr A e = and S[A] = 0, as it should be, of course, 
for any purely classical field configuration, Eq. (16 .11 . There is, however, a wider 
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class of configurations, that leave the action d6.5l > invariant (zero). Indeed, any 
field configuration of the form: 



Q = To AoT 



(6.10) 



where Ttvix) = %-t'< an d thus commutes with Qq, obviously does not change 
the action (16. 5> . This is the zero-mode Goldstone manifold. The standard way 
to introduce the massless modes ("spin-waves") is to allow the deformation ma- 
trices T to be slow functions of t + 1' and r. Thus the expression (16. lOt param- 
eterizes the soft modes manifold of the field Q. One may thus restrict oneself 
only to the field configurations given by Eq. J6. lQi and disregard all others (mas- 
sive modes). An equivalent way to characterize this manifold is by the condition 
(cf. Eq. (gij): 



Q =1. 



(6.11) 



Our goal now is to derive an action for the soft-mode field configurations 
given by Eqs. J6. lOi or ( 16.1 1> . To this end one substitutes Q = T o A o T _1 
into Eq. ( 16. 5> and cyclically permutes the T matrices under the trace sign. This 



way one arrives at T 1 o Q Q oT = Q Q +T 
[dt + vpWr, T]-, where one has linearized the dispersion relation near the Fermi 
surface k 2 /(2m) — /i w vpk — > ivpW r . As a result, the desired action has the 
form: 



iS[Q] =Trln 1 + iQT -\d u f ]_ + iQT ~ x [v F V r , t]_ 



where Q is the impurity dressed Green function, defined as: (Q Q 
For practical calculations it is convenient to write it as: 



(6.12) 



1. 



Qe{k)-- 

with 







1 



= 2 ( Qfik) [T + A e ] + G£{k) [1 - AJ ) , (6.13) 



(6.14) 



Notice that J2k^( k ) = and Efc^fW^C*) = 2?ri/r ' while the 

other combinations vanish: Y, k G?{ k )G? ( k ) = J2 k SH k )Gf(k) = 0, due to 
the complex £fc -plane integration. 

One can now expand the logarithm in Eq. (16. 12t to the first order in the <9 t term 
and to the second order in the V r term (the first order term in V r vanishes due to 
the angular integration) and evaluate traces using Eq. ( I6.14l >. For the dt term one 
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finds: 7rz/Tr{AT- 1 [<9 t ,T]_} = irv Tr{<9 t <3}, where one used that Tr{9 t A} = 0. 
For the V r term, one finds: — \-kvDTx{{\7 r Q) 2 }, where D = v^r/d is the 
diffusion constant and d is the spatial dimensionality 8 . Finally, one finds for the 
action of the soft-mode configurations [14]: 

S[Q] = mv Jdr Trji D (V r Q(r)) 2 - d t Q{r) - iV a j a Q(r) - ^V T , (6.15) 

where the trace is performed over the 2x2 Keldysh structure as well as over 
the N x N time structure. In the last expression we have restored the source 
term from Eq. ( 16.51 . The last term, V 1 <j{V is the static compressibility of the 
electron gas. It originates from the second order expansion of Eq. d6.5l > in V, 
while keeping the high energy part of the Q R Q R and Q A Q A terms. Despite of 
the simple appearance, the action ( I6.15> is highly non-linear due to the condition 
Q 2 = 1. The theory specified by Eqs. (I6.1H and (16. 151 is called the matrix non- 
linear a-model (NLerM). The name came from the theory of magnetism, where 
the unit-length vector, cr(r), represents a local (classical) spin, that may rotate 
over the sphere a 2 = 1. 

6.3. Usadel equation 

Our goal is to investigate the physical consequences of the NLctM. As a first step, 
one wants to determine the most probable (stationary) configuration, Q t (r), on 
the soft-modes manifold, Eq. ( 16. lit . To this end one parameterizes deviations 
from Q t (r) as Q = 7~ o Q o T -1 and chooses T = e w , where VV t .t' ( r ) is the 
generator of rotations. Expanding to the first order in VV, one finds: Q = Q + 
[Wo, Q]-. One may now substitute such a Q-matrix into the action (16. 151 and 

require that the term linear in W vanishes. This leads to the saddle-point equation 
for Q. For the first term in the curly brackets on the r.h.s. of Eq. d6 . 1 5 i one 
obtains: iTr{Wo V r D(V r QoQ-Qo V r Q_)} = -Tr{Wo V r D(Qo V r Q)}, 

* 2 

where one employed that V r Q oQ + Qo V r Q = 0, since Q =1. For the 
second term one finds: Tr{VV t t ' + dt)Q t , f }- It is written more compactly 
in the energy representation, where <9 t — > — ie, and thus the second term is: 
— iTr{W o [e, Q]-}. Demanding that the linear term in W vanish, one finds: 

V r (DQoV r Q)+j[e,Q]_=0. (6.16) 

8 One uses that v F = k/m and ^^( k )^Sf(k)^ = 2irvTV 2 F /d = 2wuD, while the 
corresponding Ft — R and A — A terms vanish. Employing Eq. 16.131 . one then arrives at i Tr{[l + 
A e ](f- 1 V r T)[i - AJ(t-lV r r)} = -iTrlCV^TA.T- 1 )) 2 } = -| Tr{(V,-Q) 2 }. 
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This is the Usadel equation for the stationary Q-matrix, that must also satisfy 

» 2 - 

Q = 1. In the time representation i[e, Q\- — > — {d t , Q}+- 

If one looks for a solution of the Usadel equation (16. 16t in the subspace of 
"classical" (having the causality structure) configurations, then the condition 

" 2 * 

Q =1 restricts the possible solutions to A, Eq. J6.9i (with a yet unspecified 
distribution matrix Tt,t' (?"))• Therefore, in the non-superconducting case the Us- 
adel equation is reduced to a single equation for the distribution matrix Tt.v (r). 
It contains much more information for the superconducting case (i.e. it also de- 
termines the local energy spectrum and superconducting phase). Substituting 
Eq. J6.9I into the Usadel equation j6.16> . one finds: 

V r (DVrJ 7 ) + i[e, T]- = . (6.17) 

Finally, performing the time Wigner transform, J-"t,t'{f) —> ^F(T,e;r), as ex- 
plained in section |3~51 one obtains: 

V r (DV r f F ) - d T f F = . (6.18) 

This is the kinetic equation for the fermionic distribution function fp(r, e; r) of 
the disordered system. It happens to be the diffusion equation. Notice, that it is 
the same equation for any energy e and different energies do not "talk" to each 
other (in the adiabatic case, where the WT works). This is a feature of non- 
interacting systems. In the presence of interactions, the equation acquires the 
collision integral on the r.h.s. that mixes different energies between themselves. 
It is worth mentioning that elastic scattering does not show up in the collision 
integral. It was already fully taken into account in the derivation of the Usadel 
equation and went into the diffusion term, D V 2 .. 

As an example, let us consider a disordered one-dimensional wire of length L 
[23], attached to two leads, that are kept at different voltages. There is a stationary 
current passing through the wire. We look for the space dependent distribution 
function, fp(e;r), that satisfies DV^f — in a stationary setup (for a space 
independent diffusion constant, D). As a result, 

f F (e;r) = f L (e) + (f R (e) - f L (e)) V - , (6.19) 

where //,(m(e) are the distribution functions of the left and right leads. The 
distribution function inside the wire interpolates the two distributions linearly. 
At low temperatures it looks like a two-step function, where the energy separa- 
tion between the steps is the applied voltage, eV, while the height depends on 
position. Such a distribution was measured in a beautiful experiment [23]. Com- 
paring equation J6.18t with the continuity equation, one notices that the current 
(at a given energy e) is given by J (e) = DV r f F = D (f R (e) - f L (e)) jh. And 
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thus the total current is J = eJ2 k J ( £ ) = eE T f de (fn( e ) - = e 2 ^V. 

This is the Drude conductivity: an = e 2 vD. 

6.4. Fluctuations 

Our next goal is to consider fluctuations near the stationary solution, Q t ,( r )- 

We restrict ourselves to the soft-mode fluctuations that satisfy Q 2 = 1 only, and 
neglect all massive modes that stay outside this manifold. As was already stated 
above these fluctuations of the Q-matrix may be parameterized as 

Q = e - w oQoe w . (6.20) 

The part of W that commutes with Q does not generate any fluctuations, therefore 

one restricts W to satisfy: WoQ + QoW~Q. Since Q may be diagonalized 
according to: 

z=(l -Ml -')■ 

any generator W that anticommutes with Q may be parameterized as 

w= (o -ij°U o)°{o -i) = {-w -wot J' (6 - 22) 

where Wt,t' ( r ) an d w t .t< ( r ) are arbitrary Hermitian matrices in time space. One, 
thus, understands the functional integration over Q as an integration over Hermi- 
tian w and w. The physical meaning of w is a deviation of the fermionic distribu- 
tion function, T, from its stationary value. At the same time, W has no classical 
interpretation. To a large extent it plays the role of the quantum counterpart of w, 
that appears only as the internal line in the diagrams. 

One may now expand the action, Eqs. (16. 151 . in powers of w and w. Since Q 
was chosen to be a stationary point, the expansion starts from the second order. 
In a spatially uniform case one obtains: 

iS^ [W] = 2nv JdrJJ w eie2 (r) [~DV 2 r + i(ei - e 2 )] w e2£l (r) . (6.23) 

The quadratic form is diagonalized by transforming to the momentum represen- 
tation. As a result, the propagator of small Q-matrix fluctuations is given by: 

{w e2ei (q)W e3ei (-q)) w = ~ ^f' 2 ^ = J-^^± DM , (6 .24) 

2-kv Dq z + ltd 2ttv 
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where ui= t\ — t2 and the object D(uj, q) = D(e±—€2,q) = (Dq 2 +i(ei — £2)) _1 
is called a diffuson. It is an advanced (retarded) function of its first (second) en- 
ergy argument, £1(2), (or correspondingly £1(2)). The higher order terms of the 
action's expansion describe non-linear interactions of the diffusons with vertices 
called Hikami boxes. These non-linear terms are responsible for the localization 
corrections. If the distribution function T is spatially non-uniform, there is an 
additional term in the quadratic action zS^pV] = — 2kvDTv{w\7 r Tw\/ r J-}. 
This term generates non-zero correlations of the type (ww) and is actually neces- 
sary for the convergence of the functional integral over w and w. In the spatially 
uniform case, such a convergence term is pure regularization (the situation that 
was already encountered before). 

One can now derive the linear density response to the applied scalar potential. 
According to the general expression, Eq. J5.28L the retarded response is given 
by 



n R (t,t';r,r') 



5 2 Z[V] 



2 6V cl {t';r')dV q (t;r) 



(6.25) 



where the angular brackets stand for the averaging over the action (16. 1 5I >. In the 
Fourier representation the last expression takes the form: 



W R {u J -q) = V +-{^f 



dede' 



<Tr{7«Q e , e+( ,( 9 )}Tr{7 c/ Q £ ' + ^e'(-g)}).(6.26) 



Employing Eq. (16.221 . one finds the linear in W terms: 



Tr{7«Q e , e+w (q)} ~ 2 (T e W^ +UJ (q) 
Tr{ 7 ci Q £ '+^(g)}~2(^ + . 



{q)T t i - W e >+u,,e>(q) 



(6.27) 



For a spatially uniform distribution (ww) = and only the last term of the last 
expression contributes to the correlator. The result is: 



1 



Dq 2 



Dq 



Dq 2 



) (w e < +u)j(L >(-q)We,e+u(q)) 

— , (6.28) 



where we have used the fact that for any reasonable fermionic distribution T±oo — 
±1 and therefore J dt{F t — T t +u) — —2u>. The fact that H(ui, 0) = is a conse- 
quence of the particle number conservation. One has obtained the diffusion form 
of the density-density response function. Also notice that this function is indeed 



64 



A. Kamenev 



retarded (analytic in the upper half-plane of complex of), as it should be. The 
current-current response function, K r (uj; q) may be obtained using the continu- 
ity equation qj + uig = and is K r (lu; q) — ui 2 U R (ui; q)/q 2 . As a result the 
conductivity is given by 

2 _ ■ 

a{w;q) = — K R {io;q) = e 2 vD — : . (6.29) 

ilu Dq l — ico 

In the uniform limit q — > 0, one obtains the Drude result: a(uo\ 0) = e 2 vD. 



6.5. Spectral statistics 

Consider a piece of disordered metal of size L such that L^> I, where / = vpr is 
the elastic mean free path. The spectrum of the Schrodinger equation consists of 
a discrete set of levels, e„, that may be characterized by the sample-specific DOS, 
v i e ) ~ En — e ")- This quantity fluctuates wildly and usually cannot (and 
need not) be calculated analytically. One may average it over the realizations of 
disorder to obtain a mean DOS: ^(e). The latter is a smooth function of energy 
on the scale of the Fermi energy and thus at low temperature may be taken as 
a constant v{(-f) = v. This is exactly the DOS that was used in the previous 
sections. 

One may wonder how to sense fluctuations of the sample-specific DOS v(e) 
and, in particular, a given spectrum at one energy e is correlated with itself at 
another energy ef. To answer this question one may calculate the spectral corre- 
lation function: 

R(e,e') = KeHe 7 ) - v 2 . (6.30) 

This function was calculated in the seminal paper of Altshuler and Shklovskii 
[24] in 1986. Here we derive it using the Keldysh NLtrM. 

TheDOSisKe) =iJ2M(e)-g£(e))/(tor) = «Vi&Htfe&»/(27r) = 
— i/>(73^>/(27r), where the angular brackets denote quantum (as opposed to disor- 
der) averaging and the indexes are in Keldysh space. To generate the DOS at 
any given energy one adds a source term — J de/ (2tt) J e J rfri/> e (r)o'3^) e (r) = 

— fj dtdt' J dr r ip t (r)Jt-t'd'3^t'(r) to the fermionic action. Then the DOS is ob- 
tained by v{e) = 8Z[J]/SJ e . After averaging over disorder and changing to the 
Q-matrix representation in exactly the same manner as above, the source term 
is translated to ttv J de/(2ir)J e J cir Tr{Q e , e (V)<73}. The derivation is the same 
as the derivation of Eq. ( I6.15> . It is now clear that u(e) = i^(Tr{(2 e e <73})Q. 
Substituting Q e ^ t = A e one finds u(e) = v, as it should be, of course. It is 
also easy to check that the fluctuations around A do not change the result (all the 
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fluctuation diagrams cancel due to the causality constraints). We are now in the 
position to calculate the correlation function: 



J(Tr{Q e)£ a 3 }Tr{4',e^3})Q - 1 



(6.31) 



Employing the parameterization of Eqs. ( 16 .201 —1 16 . 22i . one finds up to the second 
order in W : 

Tr{Q<T 3 } — 2[1 + J 7 ow +W o J 7 + w oW +w o w] . (6.32) 

Since (ww) = 0, the only non-vanishing terms contributing to Eq. (I6.3H are 
those with no w and w at all (they cancel v 2 term) and those of the type (wwww) . 
Collecting the latter terms one finds: 

R(e,e') (6.33) 

= V 2 Jdr IJ ^^ 2 ((We,elWel,e+W£,flW£l,£)K^ e2 W e2 , £ '+W £ ^ e2 W' e2 , £ '))Q ■ 

Finally, performing the Wick contractions according to Eq. J6.24I and taking into 
account that J dc\D 2 {e — t\\ q) — 0, due to the integration of a function that is 
analytic in the entire upper half-plane of e\, one finds: 

R(e, - [^ 2 ( £ - e '; ?) + D 2 (e' - e; q)} , (6.34) 

where the g-summation stands for a summation over the discrete modes of the 
diffusion operator DW 2 . with the zero current (zero derivative) at the boundary of 
the metal. This is the result of Altshuler and Shklovskii for the unitary symmetry 
class. Notice that the correlation function depends on the energy difference u> — 
e — e' only. 

For a small energy difference u < EThouiess = D/L 2 only the lowest ho- 
mogenous mode, q = 0, of the diffusion operator (the so called zero-mode) 
may be retained and thus: B,(ui) = —1/(2tt 2 uj 2 ). This is the universal ran- 
dom matrix result. The negative correlations mean energy levels' repulsion. No- 
tice that the correlations decay very slowly - as the inverse square of the en- 
ergy distance. One may notice that the true random matrix result Rrmt(u) — 
— (1 — cos(2irui / S)) I (2ir 2 ui 2 ), where 5 is the mean level spacing, contains also 
an oscillatory function of the energy difference. These oscillations reflect dis- 
creteness of the underlying energy spectrum. They cannot be found by the per- 
turbation theory in small fluctuations near the A "point". However, they may be 
recovered once additional stationary points (not possessing the causality struc- 
ture) are taken into account [25]. The saddle-point method and perturbation 
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theory work as long as ui > 5. Currently it is not known how to work with the 
Keldysh NLerM at u) < 5. 

In the opposite limit, oj > EThouless, the summation over modes may be 
replaced by an integration and thus R(uj) = —Cd/uj 2 ~ d / 2 , where Cd is a posi- 
tive dimensionality dependent constant. This algebraic decay of the correlations 
is reflected by many experimentally observable phenomena generally known as 
mesoscopic fluctuations. 

The purpose of these notes is to give the reader a general perspective of the 
Keldysh formalism, its structure, guiding principles, its strength and its limita- 
tions. Due to space limitations, I could not include many topics of contemporary 
research interests into this introductory course. I hope to fulfill some of the gaps 
on future occasions. 

I am indebted to V. Lebedev, Y. Gefen, A. Andreev, A. I. Larkin, M. Feigel- 
man, L. Glazman, I. Ussishkin, M. Rokni and many others for numerous discus- 
sions that shaped these notes. I am sincerely grateful to the school organizers for 
their invitation. The work was supported by the A. P. Sloan foundation and the 
NSF grant DMR-0405212. 



Appendix A. Gaussian integration 

For any complex N x N matrix Ay, where i,j = 1, . . . N, such that all its 
eigenvalues, A<j, have a positive real part: SKA; > 0, the following statement 
holds: 



^]=i/n^r^ e " 3 = detA ' (A - 1} 

-ooi =1 



where Jj is an arbitrary complex vector. To prove it, one may start from a Hermi- 
tian matrix, that is diagonalized by a unitary transformation: A = U'MJ, where 
A = diag{Aj}. The identity is then easily proven by a change of variables (with 
unit Jacobian) to Wi = UijZj. Finally, one notices that the r.h.s. of Eq. iA. U is 
an analytic function of both Ay and QAy . Therefore, one may continue them 
analytically to the complex plane to reach an arbitrary complex matrix Aij . The 
identity dA. U is thus valid as long as the integral on its l.h.s. is well defined, that 
is all the eigenvalues of Ay have a positive real part. 

The Wick theorem deals with the average value of a string z ai ... z ak z\, x ■ ■ . z\ 3k 

f N \ 

weighted with the factor cxp < — ZiAijZj >. The theorem states that this av- 
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erage is given by the sum of all possible products of pair-wise averages. For 
example, 



{Z a Zb) 



1 8 2 Z[J] 



Z[0] SJ a SJ b 

1 



.7=0 



\Za\ Z<i2 %bi 



S 4 Z[J] 



Z[0] 5J ai J a2 5J bl Ji 



ho 



J=0 



(A. 2) 



: Aii&i A» 2 &a + Aji& 2 \ 2 bi > ( A - 3) 



etc. 



The Gaussian identity for integration over real variables has the form: 

JV 



N 



Z W = / 11-7= e 13 ' =- 



(A. 4) 



where A is a symmetric complex matrix with all its eigenvalues having a positive 
real part. The proof is similar to the proof in the case of complex variables: one 
starts from a real symmetric matrix, that may be diagonalized by an orthogonal 
transformation. The identity iA. At is then easily proved by the change of vari- 
ables. Finally, one may analytically continue the r.h.s. (as long as the integral is 
well defined) from a real symmetric matrix Aij, to a complex symmetric one. 

For an integration over two sets of independent Grassmann variables, £j and 
£j, where j — 1,2, . . . , N, the Gaussian identity is valid for any invertible com- 
plex matrix A: 

Z[X,X] (A. 5) 

JV JV JV 

= //J J <,</i, e « 1 del .1 e IJ 



Here and Xj are two additional mutually independent (and independent from 
£ 3 and £j) sets of Grassmann numbers. The proof may be obtained by e.g. brute 
force expansion of the exponential factors, while noticing that only terms that 
are linear in all 2N variables ^ and £j are non-zero. The Wick theorem is 
formulated in the same manner as for the bosonic case, with the exception that 
every combination is multiplied by the parity of the corresponding permutation. 
E.g. the first term on the r.h.s. of Eq. iA. 3\ comes with a minus sign. 
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Appendix B. Single particle quantum mechanics 



The simplest many-body system of a single bosonic state (considered in Chapter 
13 is of course equivalent to a single-particle harmonic oscillator. To make this 
connection explicit, consider the Keldysh contour action Eq. ( 12. Ill with the cor- 
relator Eq. ( 12. 1 3i written in terms of the complex field 4>{t). The latter may be 
parameterized by its real and imaginary parts as: 



4>{t) 



i 



(p(t) -iu q(t)) 
(p(t)+iu q{t)) 



(B. 1) 



In terms of the real fields p(t) and q(t) the action, Eq. 



, takes the form: 



S\p, q] 



dt 



pq 



2^ 



(B. 2) 



where the full time derivatives of p , q and pq were omitted, since they con- 
tribute only to the boundary terms, not written explicitly in the continuous nota- 
tion (they have to be kept for proper regularization). Equation (1^2) i s nothing 
but the action of the quantum harmonic oscillator in the Hamiltonian form. One 
may perform the Gaussian integration over the p(t) field to obtain: 



S[q] 



dt 



1-2 W 2 

-q — - q 

2 q 2 q 



(B.3) 



This is the Feynman Lagrangian action of the harmonic oscillator, written on the 
Keldysh contour. It may be generalized for an arbitrary single particle potential 

U(q): 



S[q(t)] 



dt 



\ m 2 



U(q(t)) 



(B.4) 



One may split the q(t) field into two components: q + (t) and g_ (t), residing on 
the forward and backward branches of the contour, and then perform the Keldysh 
rotation: q± = q c i ± q q . In terms of these fields the action takes the form: 



S[q c i,q q } = dt 



-2q q 



d 2 q c i 
dt 2 



U {q c i + q q ) + U(q c i - q q ) 



(B.5) 
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where integration by parts was performed in the term q q q c i- This is the Keldysh 
form of the Feynman path integral. The omitted boundary terms provide a con- 
vergence factor of the form ~ iOq q - 

If the fluctuations of the quantum component q q (t) are regarded as small, one 
may expand the potential to the first order and find for the action: 



S[qcl,q q ] = dt 



-2q q 



d 2 g c i dU (q cl ) 
dt 2 dq ci 



H + 0(q 3 q ) 



(B. 6) 



In this limit the integration over the quantum component, q q , may be explic- 
itly performed, leading to a functional <5-function of the expression in the round 
brackets. This (5-function enforces the classical Newtonian dynamics of q c i : 



2 q c i _ dU{q c i) 



dt 2 



(B.7) 



For this reason the symmetric (over forward and backward branches) part of the 
Keldysh field is called the classical component. The quantum mechanical infor- 
mation is contained in the higher order terms in q q , omitted in Eq. ( IB. 6> . Notice, 
that for the harmonic oscillator potential the terms denoted as O(q^) are absent 
identically. The quantum (semiclassical) information resides, thus, in the con- 
vergence term, iOq 2 , as well as in the retarded regularization of the d 2 /(dt 2 ) 
operator in Eq. (IB. 6\ . 

One may generalize the single particle quantum mechanics onto a chain (or 
lattice) of harmonically coupled particles by assigning an index r to particle co- 

2 

ordinates: q r (t), and adding the spring potential energy: [q r +i (t) — qr(t)) 2 - 
Changing to spatially continuous notations: tp(t; r) = q r (t), one finds for the 
Keldysh action of the real (phonon) field: 



S[<p] 



dr Jdt 

c 



■ <P 



(V r ^) 2 - U(<p) 



(B. 8) 



where the constant v s has the meaning of the sound velocity. Finally, splitting 
the field into ((p+, (/?_) components and performing the Keldysh transformation: 
ip± = ip c i ± ip q , and integrating by parts, one obtains: 



S[tp+ ,tp-]= dr dt [2<p q (v 2 s V 2 r - d 2 ) ip c i - U{v d + <p q ) + Ufa* - ip q )] .(B. 9) 



According to the general structure of the Keldysh theory the differential operator 
on the r.h.s., ( — Of + v 2 V 2 ), should be understood as the retarded one. This 
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means it is a lower triangular matrix in the time domain. Actually, one may 
symmetrize the action by performing the integration by parts, and write it as: 

M-# + « 8 2 v?)\c + M-# + «, 2 v?)V 
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